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Abstract. We introduce a scalar weighed average ( "q-average" ) acting on 
concentric comoving domains in spherically symmetric Lemaitre— Tolman-Bondi 
(LTB) dust models. The resulting averaging formalism allows for an elegant 
coordinate independent dynamical study of the models, providing as well a 
valuable theoretical insight on the properties of scalar averaging in inhomogeneous 
spacetimes. The q— averages of those covariant scalars common to FLRW models 
(the "q-scalars") identically satisfy FLRW evolution laws and determine for 
every domain a unique FLRW background state. All curvature and kinematic 
proper tensors and their invariant contractions are expressible in terms of the 
q— scalars and their linear and quadratic "local" fluctuations, which convey the 
effects of inhomogeneity through the ratio of Weyl to Ricci curvature invariants 
and the magnitude of radial gradients. We define also non-local fluctuations 
associated with the intuitive notion of a "contrast" with respect to FLRW 
reference averaged values assigned to a whole domain or time slice. The q-averages 
of local and non-local quadratic fluctuations are directly and exactly related 
to statistical correlation moments of the density and Hubble expansion scalar. 
The evolution equations for the q-scalars and suitably defined perturbations 
completely determine the dynamics of the models without the back-reaction 
correlation terms of Buchcrt's average. A rigorous formalism of exact spherical 
non— linear perturbations can be defined over a formal FLRW background state 
associated to the q-scalars, recovering the standard results of linear perturbation 
theory in the appropriate limit. We briefly explore the possible application of this 
formalism to open theoretical issues, such as the relation between the growth of 
inhomogeneity and a definition of gravitational entropy. 
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1. Introduction. 

The spherically symmetric Lemaitre-Tolman-Bondi (LTB) dust models jTJ [2J [3l [4] 
are a very popular class of exact solutions of Einstein's equations. Since these 
models allow us to study non-linear effects in self-gravitating systems by means 
of analytic expressions or mathematically tractable numerical methods, they have 
been used in the literature in a wide variety of contexts: as models of cosmological 
inhomogeneities El El HI HI HI3 HU EH HH] (see [21 IS [13] for reviews), as preferred 
test models in the effort to explain observations without resorting to dark energy 
[31 HU [T31 [TCI HTl HH Ull I201 H21 HH] (see [23] for a review), cosmic censorship 
pMl H?5] and even in quantum gravity [25] . 

The LTB models are usually described by their standard original metric variables, 
which are well suited for most of their physical and cosmological applications, as 
can be appreciated in the abundant literature (see the book reviews [21 El H]). 
However, different variables may be better suited to examine the theoretical properties 
of the models in a coordinate independent manner, for example, the covariant 
'fluid flow' scalars in the "1+3" formalism [23 [2S [23 H3 IH 132] whose evolution 
equations and spacelike constraints are equivalent to the field equations (see [33] 
for an innovative approach to the models). Another set of alternative dynamical 
variables are the "quasi-local" scalars (to be denoted as "q-scalars" ) , defined as 
weighed proper volume averages ( "q-averages" ) of covariant fluid flow scalars on 
spherical comoving domains. The q-scalars can be, either functionals defined on 
arbitrary fixed domains, or functions ( "q-functions" ) when considering the pointwise 
dependence of the average on the varying boundary of a domain. By comparing 
q-scalars with the non-averaged scalars we obtain fluctuations and perturbations, 
which are exact, not approximated, quantities. The fluctuations and perturbations 
that can be "local" when the comparison is with q-functions in a pointwise manner, 
or "non-local" if comparing local non-averaged values with the q-average assigned to 
a whole domain. The q-functions and their corresponding local perturbations have 
been applied successfully to examine various aspects and properties of LTB models: 
to construct an initial value formulation [34], to examine inhomogeneous dark energy 
sources (quintessence and the Chaplygin gas) [35J 125] , to apply a dynamical systems 
approach [371 HE] j to examine their radial asymptotics [21!] and the evolution of radial 
profiles of covariant scalars and void formation [40] . 

The relevance of applying a formalism of q-averages for the study of LTB 
models can be appreciated from the fact that averaging has become recently an 
important open topic in current research in General Relativity. Evidently, as a frame 
dependent average acting on scalars, the q-average should be examined in reference 
and comparison to the similar averaging formalism of Buchcrt [41] . which is widely 
used in the literature (see the comprehensive reviews in [321 [43] and [44] [45] for 
an alternative covariant averaging formalism acting on proper tensors, see also |46j 
for further discussion on averaging). Since Buchert's evolution equations reveal the 
presence of "back-reaction" correlation terms that convey a significant modification 
of the dynamics, this formalism has a good potential for applications in Cosmology, 
for example, to understand the effects of non-linearity in structure formation [47 , 
in the re-interpretation of observations 48 , in the possibility of explaining cosmic 
acceleration and dark energy |49j . and even in modeling dark energy sources [50j . As 
physically well motivated toy models of cosmological inhomogeneities, LTB models 
provide an ideal scenario to apply and test Buchert's formalism [5T1I521 1531 1541 [551 15S] 
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(Buchert's average for LTB models is simply the q-average with unit weight factor). 
In particular, q-functions and their local fluctuations and perturbations have already 
been employed as useful tools to examine the existence of back-reaction and "effective" 
acceleration in generic LTB models in the context of Buchert's formalism [541 1551 [55] . 

In the present article we address various unresolved issues and extend and enhance 
previous work. In particular, besides the q-functions and local fluctuations and 
perturbations that were used in [33 EH [39l 001 HH EH [56] , we now consider the full 
range of q-scalars, which includes (besides q-functions) the q-scalars defined as strict 
average functionals and their associated non-local fluctuations and perturbations. 
While some of the background material in sections 3-5 has appeared before, we 
introduce new elements and further discussion. In particular, we provide in section 4 
a more thorough discussion of the properties of all q-scalars, emphasizing their role as 
LTB scalars that behave as "effective" FLRW scalars, as they (i) satisfy FLRW time 
dynamics, (ii) mimic FLRW expressions that are widely used in the literature (see void 
models [Il[T3[Tg[T71[TS[T3[5a[lTl[23[23]),and (iii) identify for each domain a unique 
FLRW background state through Darmois matching conditions (though an actual 
matching with a FLRW region is optional, not mandatory). We define local and non- 
local fluctuations and perturbations in section 5 and examine in section 6 their relation 
with all tensorial objects characteristic of the models (summarized in Appendix 
D). We show that all proper curvature and kinematic tensors and their invariant 
contractions are expressible in terms of q-scalars and their local linear and quadratic 
fluctuations (or perturbations). In particular, the curvature tensors (Riemann, Ricci, 
Weyl and electric Weyl) and their contractions are entirely expressible in terms of the 
Ricci scalar, its q-average and and its local fluctuations. While local and non-local 
quadratic fluctuations of the density and Hubble scalar are different objects, their q- 
average is the same for every domain (see proof in Appendix C), and as a consequence, 
the q-average of these fluctuations are averages of invariant scalar contractions and are 
equal to the statistical variance and covariance (correlation) moments of the density 
and of the Hubble expansion scalar. 

Considering that local perturbations convey the deviation from FLRW geometry 
through the ratios of Weyl to Ricci curvature invariants and anisotropic (shear) to 
isotropic expansion, and bearing in mind that these perturbations and their associated 
q-functions completely describe all proper tensors and scalar contractions, it is natural 
to expect that these q-scalars and perturbations should also yield a complete and 
self-consistent system of evolution equations that fully determine the dynamics of the 
models 35 , 36, 37, 38, 39, 40 . In section 7 we derive a similar set of complete evolution 
equations based on q-averages and non-local perturbations. 

In section 8 we show that the back-reaction terms that appear in Buchert's 
dynamical equations vanish identically in the evolution equations associated with the 
q-scalars. Hence, the utilization of the q-average (at least as far as LTB models are 
concerned) leads to different implications from those of Buchert's average: instead of 
modifying the dynamics by the extra back-reaction terms, the q-scalars provide a more 
elegant description, with deeper theoretical insight, of the existing fluid flow dynamics, 
as witnessed by the results of section 6: the direct (non-perturbative) relation between 
proper tensors and their contractions with local fluctuations and between the q— 
averages of quadratic fluctuations and statistical correlation moments of the density 
and Hubble expansion. These simple straightforward theoretically significant relations 
do not occur with Buchert's averaging (at least in its non-perturbative application to 
LTB models). 
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We showed in (35] [36] that a rigorous formalism of exact spherical perturbations 
on a FLRW background emerges from the q-functions (which dehne a formal FLRW 
background) and their local perturbations. We extend this result in section 9 to non- 
local perturbations, and in sections 10 and 11 by providing a thorough comparison 
between local and non-local perturbations. In particular, we discuss the cases when 
local and non-local fluctuations and perturbations are "confined" {i.e. localized in a 
given bounded comoving domain, with or without a matching with a FLRW region), 
and when they are "asymptotic" (the domain becomes the whole time slice). Wc 
also emphasize the fact that local perturbations describe the deviation from FLRW 
conditions through the local magnitude of radial gradients of covariant scalars, whereas 
non-local perturbations describe this deviation through the familiar and intuitive 
notion of a "contrast" between local values of non-averaged scalars and reference 
average values assigned to a whole domain (or to a whole slice in the asymptotic 
limit). As a consequence, non-local perturbations are more intuitive, as the sign of 
their amplitudes corresponds to the familiar positive/negative sign that we associate 
to over/under densities. This sign is the opposite for local perturbations: over/under 
densities are negative/positive. In section 12 we show that in the linear limit both 
perturbations (local and non-local) yield the familiar density perturbation equation 
of linear theory in the comoving gauge. 

As a quick application of the q-average formalism, we consider in section 13 the q- 
averages in the context of the entropy Leibler-Kullback functional from Information 
Theory proposed by Buchert and coworkers O [57l [53 [59]. We obtain the same 
conjecture whereby a positive entropy production from this functional follows from 
a negative statistical correlation between the fluctuations of density and the Hubble 
expansion scalar (which is the average of an invariant quadratic contraction of the 
electric Weyl and shear tensors) . We show by qualitative arguments (as a preliminary 
result) that this conjecture is fulfilled in the late time evolution of hyperbolic and 
elliptic models, but not in the early time evolution when the decaying modes are 
nonzero (non-simultaneous big bang). A summary and final comments are provided 
in section 14. The paper contains four brief appendices: Appendix A summarizes the 
standard analytic solutions of the field equations, Appendix B proves that functions 
of q-scalars are also q-scalars, Appendix C contains the rigorous proof that local 
and non-local quadratic fluctuations of the density and Hubble scalar have exactly 
the same q-average for any domain, and Appendix D provides coordinate independent 
expressions for all curvature and kinematic proper tensors characteristic of the models. 

2. The dynamics of LTB models. 

The "Lemaitre-Tolman-Bondi (LTB) models" are the well known exact solutions of 
Einstein's field equations characterized by the LTB metric [TJ[H[2]: 

D/2 



where R = R(t,r), R' = dR/dr and E = E(r) (we have set G = c = 1 and r has 
length units). For a normal geodesic comoving 4-velocity u a = Sq and a dust source 
T ab — p u a u b with rest-mass density p the field equations reduce to 




(1) 




(2 a) 
(26) 



2M' = 8irpR 2 R', 
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where R = u a R_ a = dR/dt, and M = M(r) is the "effective" mass (the quasi-local 
Misner-Sharp mass function [551 EES EI] ) • 

2.1. Covariant objects. 

LTB models are characterized by the following covariant objects associated with 
u a , T ab and the projection tensor h ab = g ab + u a u b : 

• the rest-mass energy density given by (2b): p = u a u b T ab : 

• the Hubble expansion scalar: 

H= 9 -, with: 9 = V a u a = h b a V b u a = ^ + ^, (3) 

• the Ricci scalar of the hypersurfaces 3 T[t] orthogonal to u a and marked by 
constant arbitrary t, whose induced metric is h ab — gijS l a 5l with i,j — r, ip: 

• the shear tensor (a ab = V( a W{,) — (0/3)h a b) and the electric Weyl tensor (E a f, = 
u c u d C acbd ): 

a ab = S e ab , E = e ab a ab = -t ( ^ - ^ ) , a ab a ab = 6E 2 , (5) 

£ a6 = £e afc , ^ = e^ 41 = -^ + yp, ^ ab ^ ab = 6£ 2 , (6) 

where C abc d is the Weyl tensor and the eigenvalues (E and £) of cr a& and E ab 
are given in the eigenframe e a b = h ab — 3n a n b , with n a — ^Jh rr 5 r a being a unit 
4 vector orthogonal to u a and to the orbits of SO(3). 

It is evident that the dynamics of LTB models becomes completely determined once 
the following covariant scalars 

{p, U, E, £, K}, (7) 

have been computed, which can be done through the metric functions or from suitable 
evolution equations. 

2.2. Analytic solutions vs the "fluid flow" covariant approach. 

In most work in the literature using LTB models to tackle different problems and 
applications (see for example [21 [3j IU [IT]) involves the use of the metric functions 





R! 
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H 


R! 
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M 
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Y 



that follow from the analytic solutions of the field equation (2a), which are provided 



explicitly in Appendix A (see equations (A.1)-(A.3|). These solutions determine the 



dynamics of the models, and in particular allow us to compute (implicitly) quantities 
like R' and 7?', which are needed to compute the covariant scalars (171) from their forms 
in @ and @-((6]). 

However, the dynamics of the models can also be fully determined by finding the 
scalars ^ through suitable evolution equations, such as those that follow from the 
"1+3" formalism of Ehlers, Ellis, Bruni, Dunsby and van Ellst [HI CM 123 EH 133 E] : 

4-7T 

H= - H 2 - — p-2E 2 , (8a) 
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p = — 3pH, 
t = -2UT, 



£ 2 -£, 



£ = -47TPE-35 (H + T,), 
together with the spacelike and Hamiltonian constraints 

'47T 



3£§=0, 



P 



8tt 



P-/C + S 2 



(86) 
(8c) 



(9) 



(10) 



Solving this system of partial differential equations is equivalent to working out the 



models through the field equations (2a|-(26| and the solutions ( A.f |-(A.3). However, 
(|7| is not the only nor the most convenient 'fluid flow' scalar representation to handle 
the models. 



3. The q— scalars. 

As shown in 34, 35, 36j[39l|40], an alternative coordinate independent "fluid flow" 
representation for the LTB models is provided by the "quasi-local" or q-scalars and 
their fluctuations and perturbations. We provide below a rigorous definition of these 
variables. ||] 



3.1. The q-averages. 

Consider an arbitrary complete [|] time slice 3 T[t] of an LTB model with metric (JlJ 
whose proper volume element is 

dV p = ^/det{h ab ) d 3 x = R ~ R '* mid drdtidtp, (11) 



1/2 

(12) 

Every 3 T[t] can be foliated by the class of compact concentric spherical comoving 
domains diffeomorphic to X>[r&] = x[ r b] x S 2 C 3 T[t], where S 2 is the unit 2-sphere 
and the closed real interval x[ r b] = {r | < f < r^} is a segment of a radial ray with 
f = marking the worldline of the symmetry center (fixed point of the rotation group 
SO (3)). Following [SS], we define 

| Sec 35 36 for a discussion on the relation between the definition of q— scalars (q-averages and 
q— functions in ( |13| l and j!8) ) to the quasi-local or Misner— Sharp mass when A = p. The notation and 
units here are different from those of previous articles: q— scalars were denoted as averages "(A)" in 
[34| and as A* in |35II36| , while the symbol m was used for (4/3)irp and k for K in 55, 54, 56, 39, 40 . 
The notation used in | |13| l expresses the domain dependence of the q-average by the symbol "[!•{,]" that 
univocally identifies any domain T>[rb] by its boundary: the comoving 2-sphere marked by constant 
arbitrary r = rj,. We use the subindex q to distinguish the q-average from the "standard" proper 
volume average functional u (A) p [r^Y of Buchert's formalism, whose weight factor J 7 is a constant. 
Both averages coincide only for parabolic models | |A.1| for which T = 1, E = and self-similar LTB 
models for which T ^ 1, E = const, (see [56]). 

§ See [401 156] for the case of non-complete slices intersecting a curvature singularity. 



where (from (2a 



R 2 



2M 
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The q— average. Let X(D[r b ]) denote the set of all smooth integrable LTB 
scalar functions in a domain T>[r b ] with r b fixed but arbitrary marking its 
boundary. The q-average of a scalar A € X(D[r b ]) is the linear functional 
( } q [r b ] : X(V { r b ]) -> R such that: 

fi b A(f)R 2 (f)R'(f) df 



A i-> (A} q [r b ] = 



J v[rb] ATdV p 



J? R 2 (f)R'(f) df 



(13) 



where the "weight factor" T is given by (12). 



It is straightforward to evaluate in an arbitrary fixed domain T>[r b ] the q-average of 
those covariant scalars A = p, R, K, that are common with a FLRW cosmology: 



3 {P) q [n\- 3 jr bR2R , d - 



2M b 
R 3 



(n) q [r b ] 



;; b n r 2 R'df 

J rb R 2 R'df 
, K r , Q ICR 2 R'df 



R b 
R b 



Eb 

Rl 



Rl 



(14) 
(15) 
(16) 



where R b = R(t,r b ), M b = M(r b ), E b = E(r b ) (subindex b will denote evaluation at 
r = r b ), and we simply applied the definition ( 13 1 to ( p&| , Q and Q. The q-averages 
( 14 1— ( 16 ) transform the field equation ( |2a| for every T>[ro] into 

3tt 



(17) 



/qL ' u ' 3 

which is formally identical with the Friedman equation for a dust FLRW model given 
in terms of the q-averages of the equivalent LTB scalars. 

It is important to remark that the weight factor ( 12 ) is an invariant quantity, 
since R and M are invariant scalars in spherically symmetric spacetimes jSHj- In fact, 
in the Newtonian limit (2M/R <C 1) this factor reduces to the total binding energy 
of comoving dust layers, while in the special relativity limit it becomes the "gamma 
factor" EL]. 



3.2. The q-functions. 



Since r = r b in (131 and (14)-(16) is arbitrary, the definition of the q-average 
functional leads in a natural way to local functions with this corresponding rule. We 
define then: 

The q— functions. For every 3 T[i] foliated by concentric domains T>[r] with 
r > 0, the q-function associated with any A E X(T>[r]) is the real valued 
function A q : 3 T[t] — > K such that 

A q {r) = {A) q [r], (18) 

holds for every r > 0. Strictly speaking, the q-functions are q-averages that 
depend on a varying domain boundary. 

Evidently, (14)-(16) and (18) imply that the q-functions associated with A = p, R, K, 

are 



in 



Pq 



n, 



8tt 



2M 

W 



R - R 



I — J 7 

R 2 



K a 



(19) 
(20) 
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where all functions above (save those inside the integral sign) are real valued functions 
that depend on the upper bound of the integration range < r < r for arbitrary r (or 
the domain boundary for arbitrary domains). 

3.3. Functionals vs functions. 

The difference between the average functional (A) g [r;,] and its associated q-function 
A q (r) is subtle but important (see figures la and lb of references (Ml [551 [55] ) . For 
every fixed arbitrary domain T>[r b ] C 3 T[t] they are identical at the boundary r = r bl 



but they differ for all points r ^ r b because (A) g [rf,] has been assigned by (13) (as a 



functional) to the whole of and must be considered as a constant for points r < r b 
inside this domain in each time slice 7~[i], while A q (r) (as a function) varies smoothly 
along these inner points (and also along outer points r > r b ). As a consequence, the 



(A) g [rfc] and A q {r) satisfy the same local derivation rules ( 23 )— ( 24 ) at the boundary of 
each generic domain (see next section), but behave differently when integrated along 
any fixed domain in the range < f < r b . Hence: 

{{A) q [r h ]) q [r h ] = {A) q [r h ] but ( A q (r) ) q [r b ] ? A qb = A q (r b ), (21) 
holds for every scalar and every domain. The difference between functionals and 
functions is clearly illustrated in figure 2. 

34. Notation. 

We will adopt henceforth the following conventions in order to simplify the notation: 

• The time dependence of functionals and functions will be omitted unless it is 
needed to avoid ambiguities. The symbol f will be used as a dummy variable 
inside integral signs. 

• Unless stated otherwise, radial dependence given in terms of r (i.e. u M(r)" or 
"p q (r) n or u {p) q [r]") or altogether omitted (i.e. "M" or "p g " or "(p) 9 ") will 
be understood as local point-wise dependence like q-functions, whereas radial 
dependence in terms of r b will denote functionals that correspond to a fixed 
arbitrary domain T>[r b ] and functions that equate to them (as in the right hand 



sides of 1 14 H 16 1) 



4. Properties of the q scalars. 

4-1. Time and radial derivatives. 

It is straightforward to show that the following commutating rule holds for every scalar 
A 

l (A>9_ (§?) =(4,-(4=3(^),-3W,(i)„ (22) 

where we used ((aHJ) and expressed ^ as 3H = 9 = (R 2 R')/(R 2 R r ). The following 
useful relations are readily obtained: 

W' q = ?r{A) q = ^[A-(A) q ], (23) 



dr" ,q R 
1 

RW)Jo 



( A )« = A (r) ~ i / A '(f) R 3 (r) df, (24) 



where (24 1 follows by integrating (13) by parts. The same properties hold identically 
for the q-functions A q . 
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As we prove in Appendix B (see also Appendix C of all scalars expressible as 

functions of q-scalars are themselves q-scalars. As an example, consider the Omega 
factor (analogous to a FLRW Omega factor) defined as a q-function il q = U(p qi T-L q ) 
or fl q = U(p q ,]Cq) by: 



a, 



871-pq 



iir/3)p q 



(8ir/3)p q - JC q 



i 



JCq 

H 2 
' L q 



JCn 



(87r/3)/9 g - JC q 



whose associated local scalar is given by ( |B.3[ ) as 

, p~p q 2(H-H q ) 



Pq 



rL a 



(25) 



(26) 



where we used the quotation marks to emphasize that ft ^ 8np/(3'H 2 ) = U(p,'H). 
Notice that ffl can be expressed as ft = Q q (l + <5 (n) ) with 5^ given by 



4-3. The q-scalars satisfy FLRW evolution laws and mimic FLRW expressions. 

It is straightforward to show that the q-scalars (whether evaluated as q-functions or 
as functionals in fixed domains Z^t,]) satisfy FLRW evolution laws: [jj 

Pq = — SHqPq, JC q = — 2'H q }Q q 

nl - —p q , a q = n q (<n q - 1) , 



T-Lq — 



(27a) 

"q g pqi ""<? — ""<? V""9 -V t (276) 

which are identical to the evolution equations satisfied by the equivalent covariant 
scalars p, H, JC, Ct of a FLRW dust model (a tilde will denote henceforth FLRW 
objects). This fact evidently singles out the q-scalars as LTB scalars that behave 
as FLRW scalars (in the sense that they comply with FLRW dynamics). 

By parametrizing the metric functions in in terms of their fiducial values at an 
arbitrary time slice t — t and making the coordinate choice R Q — r, we can transform 
into the "FLRW look alike" metric: 

Tdr 2 

— ^ + r 2 (dtf 2 + sin 2 d&ip 2 ) , (28a) 



ds 2 = -At 2 + a 2 {t,r) 



1 



R 
Rq 



R 



F = 



R'/R 

Rq/Rq 



rR! 



= 1 



ra 
a 



(286) 



and as a consequence, the q-scalars (19), (20) and (25) take the following "FLRW 
look alike" forms: 



Pq 



PqO 
„3 ■' 



JCn 



tc 



tt — ttqO 



n 



q0 



a 

1 - Q 



q0 



n 



q0 



„0 



+ (1 - ftqo) a 1 



(29a) 
(296) 



where the subindex o indicates evaluation at t = to. The forms (28a)-(286) and 



( 29a) -(296) have been used in most recent articles looking at LTB void models 
|Ti n^ TT^ n T7tP^fT ^ [ ^ [ ^ [ ^ [^. where they are loosely introduced as "convenient" 
ansatzes without any justification other than their "FLRW look alike" forms, thus 
ignoring the fact that they can be defined rigorously as q-scalars that emerge from 



the weighed average (13). 

|| For the functionals (A) q rJ the derivatives involved are (A)' q [ rhL which can be evaluated (at r = rj,) 
either directly from ||14| H 16p , or with the commutation rule i ]22[ i and the forms of the local (non- 
averaged) scalars in \2b\ , (3L Q a nd ||26| . If using \22\ for computing ("H^lrf,] and (f2)'q[r6] we also 
need to use the identities 4lT and (|42|| that are proved in Appendix C. 
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4-4- The q-scalars associate to each domain a unique FLRW "background state". 



For every fluid flow FLRW scalar A = p, H, K,, ft and its "FLRW equivalent" LTB 
q-scalar A q = p q , H q , fC q , Q q , the value A qb — A q (t, r b ) identifies, for each domain 
T>[r b ] of an LTB model, a specific FLRW dust model that can be smoothly matched 
(under Darmois matching conditions) at an arbitrary finite comoving radius r = r b 
that also marks the boundary of 2?[r&]. Consider a dust FLRW universe with metric 



ds 2 



-dt 2 + a 2 {t) 



dr 2 



1 



(dtf 2 



'■Dd(p 2 



(30) 



where fco = ± 1 , £o is an arbitrary length scale and a(t) is the dimensionless FLRW 
scale factor. Darmois conditions (necessary and sufficient) for the smoothness of the 
matching of ( 30 ) with an LTB model along a comoving boundary r = r b are given by 

mum 



Pqb = (p)<?M = Pit) 



Po 



H q b = WqVb 

R b = a(t) r or a b 



H(t) = 

a 



ICqb = (IC) q [r b ] = K{t) = 



/Co 
la' 



(31a) 

(316) 
(31c) 

where the subindex denotes evaluation at a fiducial hypersurface t — t a and a(t ) = 1 
holds. We note that conditions (31a) and (31&) require the continuity of the A q 



8np 
W 2 



if we would consider a (hypothetical) smooth matching of a domain ~D\rb\ of an 
LTB model with a F L RW region . Th i s fea ture is strikingly evident if we use the 



parametrization (|28a|-(|286|) and (|29a 
A 



-(2961 with A = A qb0 = A q (t ,r b ) holds for 



p, H, /C, Q. However, from (35) and (57a)-(576), it is evident that we have 



A b ^ A qb and A' q , A' q ^ in general, and thus the local scalars \l\ and the gradients 



A' and A' q do not comply with the matching conditions (31aH31 

It is important to remark that the continuity of the A q under Darmois matching 
conditions is simply a formal rigorous procedure to identify for every D[r b ] of an LTB 
model a particular FLRW dust model that can be defined as a reference "background 
state". We use the term "state" to emphasize that this identification does not 
force us to consider an actual matching with a FLRW region, which would yield a 
"Swiss cheese" configuration in which the background state becomes also an actual 
background spacetime. Likewise, the A q also allow us to define a "FLRW equivalent" 



region to every domain 2?[r&], as they provide through (31a)-(31c) the values of 



the FLRW scalars A if the whole domain was replaced by an equivalent spherical 
comoving section of a FLRW spacetime. Therefore, having determined the FLRW 
background state, the discontinuity of the local scalars and the gradients A' and A' q is 
not problematic if we do not wish to construct an actual Swiss cheese model through 
a smooth matching at the domain's boundary r — r b . In this latter case we can avoid 



discontinuities by demanding (besides (31a)-(31c)) also the continuity of local scalars 
at r = r b through the following extra supplementary condition: 



A' b = 0, =*> A b = A qb for A = p, K, H, ft, 



(32) 



1 q b — u ) —f si-b — rt-qb 

which forces A to coincide with A q at r = r b , and thus explains the appearance (see 
panels (b) and (d) of figures 1 and 2) of "humps" (if A' > 0) or "bags" (if A' < 0) in 
the radial profiles of local scalars (this has been noted in the local density profiles in 
Swiss cheese models in LTB void models). 
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5. Fluctuations and perturbations of q scalars. 

5.1. Local fluctuations and perturbations. 

If A and A q = (A) q are both evaluated as real valued functions on the same arbitrary 
value r that denotes a varying boundary of concentric domains T>[r] for r > 0, then a 
local fluctuation can be defined at each r by the simple pointwise comparison: 

D(A)(r) = A{r) - A q (r) = A(r) - (A) q [r]. (33) 

By normalizing D(A) with A q we obtain the useful dimensionless relative local 
fluctuations that will be called local "perturbations": 



which comply (from ([23]) and (24 1) with the following useful relation with radial 



gradients of A q and A (also valid for the (A) qJ 



5 {A) = ^ = i TzkzTT, I A\f) R*(f) df, i 3* , 



ZR'/R A q (r)R3(r)J 



that lead, using (23), (24) and (B.3), to the following linear algebraic relations among 
the flW; 

25™ =n g S^ + [l-fi,]#*>, (36) 
tfW = j</0 - 26™ = (1 - n 9 ) - 6™) , (37) 



where f2 9 is given by (25 1 and 5^ above is consistent with Q = O g (l + (5^-*) in (26) 



5.2. Non-local fluctuations and perturbations. 

As opposed to local fluctuations that compare A with A q = (A) q at the same r, we 
can define for every fixed domain D[r b ] non-local fluctuations 

D NL (A)(r, r b ) = A(r) - (A) g [r b ], (38) 

that compare local values A(r) inside the domain with the q-average (functional) of 
A, which is a non-local quantity assigned to the whole domain (notice that at every 
3 T[t] the value (A) q [r b ] is effectively a constant for all r < r b and a function of t for 
varying 3 T[£]). Non-local perturbations are readily defined by 

^\r,n) = ^^ = ^§^ , 0<r<r b , { A) q[ r b] ^0, (39) 



and, evidently, do not comply with ( |35[ ) and the properties that follow thereof (notice 
that d/dr[6if(r,r b )} = A' / (A) q [r b ]) . 

5.3. Statistical fluctuations. 

Since evaluating a q-average for a fixed arbitrary D[r b ] involves integration through 
the range < f < r b , and the functional (A) g [rfc] is effectively constant in this 



range at every 3 T[i], then, following (21 ), the q-average of local and non-local linear 
fluctuations is different: 

(B(A)) q [r b }^0 but (D NL (A)),[r fc ] =0, (40) 

which means that the D NL (A) are statistical fluctuations, but the D(A) are not. 
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However, the average of quadratic combinations of local fluctuations of p and % 
coincides with the average of the same combinations of non-local fluctuations. As 
we prove in Appendix C (see also Appendix B of 56 for similar results in Buchert's 
averages), the following results are valid for A, B = p, % and for any domain Z>[r&]: 



XA) D NL (E) ), = (AB) q - (A) g (B} q = Cov 9 (A, B), 



(41) 
(42) 



( [T>(A)Y ) q = ( \D m (A)Y ) q = (A 2 ) q (Ay q = Var g (A), 
(D(4)D( J B)> g = (D NM 
where we omitted the [r b ] symbol to simplify notation (A, A q inside inside the ( ) q 
brackets depend on f < rb), whereas Var g (A) and Cav q (A) are the quadratic variance 
and covariance (correlation) moments associated with (A) g [rb]. 



6. Relation of q scalars and their fluctuations with curvature and 
kinematic invariants and tensors. 



As shown in Appendix D, all proper curvature and kinematic tensors characteristic 
of LTB models are expressible in terms of algebraic combinations of g ab , u a and the 
eigenframe e ab defined below ([5| and ([6]), with their scalar coefficients given entirely by 
four scalar invariants: the Ricci scalar (TV), the Newman-Penrose conformal invariant 
($ 2 = —£), the Hubble expansion % — 9/3 and the eigenvalue of the shear tensor, £ 
(see ^ and ||5j)). As a consequence, the q-scalars p q , H q are coordinate independent 
quantities, as they are directly related to these four scalar invariants: 
4-7T 1 

—p q = -n-y 2 , n q = n + j:, (43) 

where we used @, (|3j)— @, ([l4j)-([l7]) and/or ([l9])-(|20j), and the fact that K = 8irp. 
We remark that these relations hold for their associated q-averages (p) q , {T~L) qi {K<)q 
as functions of a varying domain boundary r. The expressions for other q-scalars, JC q 

S, W can be found by eliminating these q-scalars from the 



and Q q , in terms of 1Z 



constraints (20) and (25) and substituting (43). 



It is easy to show that the eigenvalues of the shear and electric Weyl tensors are 
local linear fluctuations over the q-averages (or q-functions) of W and p: 

E = --D(H) = -[H-H q ] £ = -^D(p) = -y[p-p 9 ], (44) 



where we applied (19) to ([3]), J5| and Considering that *2 = —£ and using (D.4), 



the curvature and kinematic proper tensors (see Appendix D) are expressible in terms 
of p, H and their fluctuations D(p), D(H): 



ab 8n fnxl a X b ] i axl- 
ed = -J p \ 1° 4 I' 



< b] 



°[c°d] 



47T 

T 



D(p) ( 



hi" - 3mM a ) e b J, 



K a b =4np(h a b 

47T 



vrub) 



D(p) e afc 
o 

a ab = -B(H)e ab , 



s~iab 



~T>(p) (/ 



T~L a b 



47T. 

Hh ab --D(H)e ab , 



while their scalar contractions take the form: 

2567T 2 ( 5 9 N 
[D(p)] 2 ' - 



abed 



CabcdC 



abed 



3 

256tt 2 



P(p)] 



a ab a ab = 6[D(H)] 



O-abE 



4 P , 

Q 171 J7&b 

oHj ab tj 

4tt 



3ur r U 1 ' 



K ab K ab = M^p 



2 „2 



ab 



D(p)D(H), 



(45 a) 
(456) 
(45 c) 
(45 d) 

(46 a) 

(466) 
(46 c) 
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which depend exclusively on p, H and their fluctuations. Con sider in g tha t we have 
A = A q (l + 8^) and D(A) = A q 5 {A) for all A, the tensors ( |45a[ )^(45d| and their 
contractions ( |46a[ )-(46c) can be entirely given in terms of p q: H q and their local 
perturbations: 

Air 

E a b =--p q 8^e a bl 



Cub 
^,1 



77 a b 
IK -cd 



477 

T 



8,7 



n a b =4mp q {l + 8^){h a b +u a u b ), 



r a b =-H q 5Wel 



abed 



T^-abcdK 



Cr^abcd 
abed^ 



2567T 2 

3 

256tt 2 



nja nj 

it,, — n„ 



(P)l2 



p\ ft 



h" 



\1 



[Pq 



0~abV 



ab 



4 

o rp rpab 

nj njab nj'. 

TtabTt — >l ( 



l + (2 + 3*< w >)*M 



(47a) 
(476) 
e fc] 

e d]> 

(47c) 
(47d) 

(47e) 

(47/5 



which shows that the variables {p q , Hq, S^ p \ 5^} should provide a full dynamical 
representation for the models, so that 5^ — 5^ = indicates a FLRW limit 
characterized by the vanishing of the shear tensor and Weyl part of the curvature 
tensors. Notice that all these tensors can also be expressed in terms of other q-scalars 
(/Cq, f2 q ) and their local perturbations (6^, 5^), as the latter can be obtained from 
p q , Hq, 5( p \ 5^ by means of the algebraic constraints (20), (25), (36) and (37). 

Bearing in mind that 7Z q = &Trp q , we can express the conformal invariant *2 as a 
local fluctuation of the Ricci scalar: *2 = (1/6)D(7*L) = (47r/3) D(p), and thus express 
curvature tensors (the electric Weyl, Weyl and Riemann) and their contractions in 
terms of the the Ricci scalar and its fluctuations: 



E ab = — -D(7£) e ab , 
o 



(-tab 



K 



3 V I c 



T^-abcdK 



abed 



«,/] 
5 



*D(7e) (hi - 3u [c u^ 



d lc°d] 



TV 



6 

CabcdC 



abed 



[T>{TL)f 



(48 a) 
(486) 
8E ab E ab , (48c) 



where we used (48a) and (D.3). These expressions provide an elegant geometric 
interpretation of the inhomogeneity associated with the Weyl part of curvature tensors 
and scalars in terms of fluctuations of the Ricci scalar. 

Since a non-vanishing shear and Weyl curvature are indicative of inhomogeneity, 
it is very useful to express the basic local perturbations 5^ and 5^ as 

£ 6* 2 £ 



P~ Pq 
Pq 

ft a 



S (P) 



5 (H) _H-H q 



iTVPg/3 

£ 

T~Ln 



n-m 2 
s 

•H + S " 



C 

i + C 



where the quotients 



6*2 

K ' 



(49 a) 
(496) 

(50) 
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provide an invariant measure of inhomogeneity through the ratio of Weyl to Ricci 
scalar curvatures (£) and the ratio of anisotropic to isotropic expansion (£) |^] . The 
remaining pertu rba tions 5 ^ and 5^ can be expressed in terms of £ and £ by using 
the constraints (36 1 and (37) to eliminate them in terms of 8^ and 6^ and then 



substituting (49a) and (496) 



Since E and £ are local linear fluctuations of rl and p, we have for all tq: 



(s),W^o, 



(51) 



However, as a consequence of (41), (42) and (44), the averages of the scalars E , £ 



and E £ are averages of local quadratic fluctuations, leading to the following statistical 
fluctuations that relate to the variance moment with respect to (rl) q and {p) q at any 
domain: 



(E 2 ) 9 = ( \D(H)] 2 ) q = {H 2 ) q - {U)l = Var 9 (H), 



<£E) 



47T 

T 

4tt 



-in 

T 

47T 



[(A - (P) 2 q ] 



4-7T 

T 



Var 9 (p), 



47T 



[( P H) 9 -(p) 9 (H) 9 ] = — Cov 9 (p,H), 



(52 a) 
(526) 

(52c) 



so that, in view of ([5| and (46c), these fluctuations relate to averages of quadratic 
contractions of the shear and electric Weyl tensors: 



{<J a bO- ah ) q 



(E a bE 



ab\ 



= 6(£ 2 ) 9 



6Var„(H), 

= 6((* 2 ) 2 ) g 



327i- 



Var g (p), 



(<JabE ab ) q = 6(E£) 9 = 8nCav q (p,H), 

while the q-averages of contractions of the Riemann and Weyl tensors are expressible 
either in terms of p or 1Z and their statistical variance moments: 



(53 a) 
(536) 
(53c) 



(R-abcdTl 



abcd\ 



(CabcdC 



abcd\ 



256tt 2 
3 

256tt 2 



Var,(p) 
Var g (p) = 



5 

-,</'> + 5 
4 



(P 2 ) q 



Var g (ft) 



Var g (ft) 

{E a bE ab ) q 



-(K 2 



(54a) 



(546) 



where we used the fact that (TZ) q = 8n{p) q and !Z a blZ = 1Z . We can express the 
quadratic ratio of Weyl to Ricci curvatures as a sort of "standard deviation" of p with 
respect to (p) q : 



6(E ab E ab ) q _ 6((* 2 ) 2 ), _ 6Var 9 (p) _ (p 2 ) q - {p) 2 q 



(pi\ 



(p*\ 



(55) 



(TZ a bTZ ab ) q ((R.) 2 )q \h> lq \y Iq 

A similar standard deviation of H with respect to (H) q follows as the quotient of 
averages of quadratic covariant scalars o- a bO~ ab and H 2 — 9 2 /9: 



(c- ab c- ab )g Var g (H) = (H 2 ) q - {U)\ 

6(w) q (w) q m q ■ 

where we used |5]) and (44). 



(56) 



1 Wainwight and Andrews (reference [33]) also emphasize the role of the ratio of Weyl to Ricci 
curvature as an invariant measure of inhomogeneity in LTB models. However, their study was not 
based on averaging and they consider a quadratic positive definite ratio C a bcdC cd / r R- a bT^- ab ■ As a 
contrast, the sign of first order ratio ^2/^ that we use here is not a priori defined. 
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7. Evolution equations. 



While the q-scalars A q or (A) g [r&] satisfy FLRW evolution laws, such as (27a)-(276), 
these scalars are not fully determined by these FLRW equations because, unlike their 
equivalent FLRW scalars A, they either depend directly on r or r&. The missing 
dynamical information is provided by the evolution equations for the 5^ and the 

§ (A) 



7.1. Local perturbations. 

Since all the local covariant scalars ^ (as well the tensors and scalar invariants in 
(47a|-(47<7)) can be expressed in terms of q-scalars and their local perturbations: 



Pi 



1 + 



rt — l~Ln 



i + 5P0 



fC = /C 



1 + 



£ 



4tt 



Pq 



5^, 



(57a) 
(576) 



the evolution equations for the variables p q , H q , S^ p \ 5^' should yield a self- 
consistent and complete set of evolution equations for the LTB models. This system 
follows readily by inserting (57a I and |576| into the "1+3" system (8a|-(8d) and its 
constraints (|9])-(10). The result is the evolution equations 

(58a) 



Pq 



— 3p q Hq, 
n,2 47T 

— n q 
-3(1 



'Pq, 



m = (1 + 3^)) H g 5W + %^(5M ~ 

Sri a 



plus the algebraic constraints 
8tt 



172 

H q — 



~Pq - & 



28M = n q 6<p) + [l - n,] <5 (fc) , <5 (0) - 5^ - 26M , 



(586) 



(58 d) 



(59) 



that exactly coincide with the general relations (20), (36) and (37), hence they hold 
at all t [i.e. they propagate in time). The following points are worth remarking: 

• The constraints Q of the 1+3 system (8a)-(8d) are satisfied trivially from (35) 
applied to <5 (p) and 5^. 

• The hrst two constraints in ( 59 ) follow from substituting ( 576 1 into ( 10 1 and using 
(2a) and (19), while the third one is obtained by differentiating (25) with respect 
to r and applying (35 ). Once (58a)-(58d) is solved these constrains and (25 1 allow 
us to compute Q q , K q , 5^ k \ S^ l> . 

• The fact that the constraints (59) of the system (58a)-(58d) are algebraic implies 
a great simplification of the numeric treatment of these fluid flow evolution 
equations, as they can be effectively integrated as a system of autonomous ODE's 
in which the initial conditions are restricted by the algebraic constraints, and thus 
it is far easier to handle than the 1+3 system (8a)-(8d) in which the constraints 
([9]) are partial differential equations on r that must be solved before dealing with 
the time equations. 
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Considering the relation between the A q and 8^ in the constraints (59), each of the 
following combination of four scalars: 



any two of A q = p q , H q , JC q , Cl q , 
plus 

any two of 5™ = 5™, S™, 5™ , 8™, 



(60 a) 
(606) 



provides a full covariant scalar representation for the models, since the remaining 
pairs of A q and 8^ can be obtained from these algebraic constraints. Each of these 
representations yields a self-consistent and complete set of evolution equations that is 
alternative (and equivalent) to that of local covariant scalars given by ([7|, and thus, 
completely determine the dynamics of the models. 

The evolution equations (58a|— (58d| correspond to the representation 



{p q , H q , 5 (p \ 6W}, which is useful to compare with the spherical collapse model 
and perturbative scenarios of structure formation that consider the density and Hub- 
ble velocity as dynamical variables. However, a more appropriate representation 
for cosmological applications (for example void models) is furnished by the scalars 
U» 6M, 6M}. 



7.2. Non-local perturbations. 



If we consider non-local perturbations ( 39 1 in an arbitrary hxed domain T>[r b ] , then the 

only local scalars in ([tJ) that can be expressed in terms of the variables {(^4) 9 [rh], 6^} 
in a similar manner as in J 5 7 all are: 



p = (p) q [n\ 



1 



s (p) 



n = (H) q [r b ] 



1 



K. = (IC) q [r b ] 



1 + S {K) 

1 T "NL 



(61) 



expressed as (576 1 



The remaining local scalars, £ and £ , given by local fluctuations by (44), cannot be 

4tt 



E(r) = H q (r) - (H) q [r b ](l + 5™), £(r) = - p q {r) - (p) q [r b }8 



dp) 



(62) 



and as a consequence, the variables {(j4) g [rh], 8^}, A = p, %, K. do not provide a 
complete scalar representation of the dynamics of LTB models, which is not surprising 
because the (-A) 9 [rf,] depend only on t and constitute boundary conditions for the A q 
at r = r b . From (39) and (61), local and non-local perturbations are related by 



1 + ^ = 7^(1 + ^), 
( A )g[n] 



(63) 



which shows that they only coincide at the boundary of each V[r b ] where A q (r b ) = 
(A) g [rt]. Since their time derivatives are evidently different, it is reasonable to expect 
that the evolution equations of the 8^ and 5$ wu l be different. Inserting (63) for 



A = pfH into (58c|-(58d| yields these equations 



3(l + ^ p L ) )(H) 9 [r & ]^ ) , 



8& = - (1 + 38 



^)(H) q [r b ]8^ + ^^(8^-8^ 



-2{H) q [n][l- 



(n) q [r b ] 



+ MH q (r)-(H) q [r b })8^\ 



(64a) 



(646) 
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which, in order to render a fully complete system to describe the dynamics of the 
models in an arbitrary fixed domain Z>[r&], needs to be supplemented by the evolution 
equations for (p) q [r b ], (H) q [n] and p q) % q : 

4-7T 



(p)'q[n] ■ 



-{U)t[r b ] 



-{p)q[n] 



-HH) q [n}(p) q [ro} 



Pq( 



U q (r) = 

) = -m 



-K{r) 
'r)p q (r) 



4tt 



Pq{r), 



(65 a) 
(656) 



while the constraints take the form (59) and with the expressed in terms of the 



by (63 1 (there is an extra constraint given by (jlT])). It is interesting to remark 



that (64a) is identical to (58c), but (646) differs from rt58c|) by the terms with % 



<~qi 

which explains the need to add the evolution equations for p q and H q . Notice that the 
non-local evolution equations become identical to the local ones (58a)-(58d) at the 
domain boundary in the limit r — ¥ r b where H q (r b ) = (H) q [r b ]. Also the non-local 
evolution equations depend on averages p q (r) = (p) q [r], H q (r) = (H) q {r] for inner 
points of V[r b ] and are considerably more complicated than (|58a|)-(58d). The self- 



consistency of the system (|64a 
derivatives ]' 



(656) can be proved easily by comparing the mixed 
obtained from = A' (r)/(A)^ 



in 



with those that follow from 



the radial derivative of the right hand sides of (64a)-(646) 



8. Comparison with Buchert's averaging. 



The most important difference between the q-average ( 13 ) and the standard proper 
volume average of Buchert's formalism (the average (13) with J- = 1) is the 
fact that the evolution equations that follow from ( 13 ) and its fluctuations lack 
the "back-reaction" correlation terms that appear in Buchert's evolution equations 
HUE! S3]. This fact follows readily by comparing ([65a|-([656| (or (|58a|)-(([58b|)) 
with the equivalent equations in Buchert's formalism. 



Using (44) to express £ in terms of H — T-L q = D(?l) and the commutation rule 
(H)' q — (fi.) q taken from (22), the weighed average (13) applied to the Raychaudhuri 
equation 



yields: 



Anr 

(H) q [r b \ = -{H) 2 q [r b } - —{p) q [r b ]+2Q q[n _ 



(66) 



where Q q [r b ] is the back-reaction term 
Q q [r b ]^{{H-{H) q f) q -{{H 



n q f), 



[d nl (w); 



[rb 



in 



o. 



(67) 



which vanishes identically for every domain as a consequence of (41 ) applied to A = % 



(see ( C.l )-( C.3 ) for the proof of this result, see also [56]). Hence, ( 66 ) reduces exactly 
to the FLRW Raychaudhuri equation (65a) with zero back-reaction. 

It is straightforward to prove that expressing the 1+3 local scalars ^ as ( 61 ) 
and (62) and applying the time derivative (22) to the Raychaudhury equation (8a) 



and to the remaining 1+3 equations (86)-(8d), yields, for any domain V[r b ], the full 
set of evolution equations (64a)-(656) without any back-reaction term (just as the 
substitution of ( |57a[ ) and (576) into (86)-(8d) yields the system (58a|-(( (58d | without 
back-reaction terms). Hence, the weighed q-average (13), either with local or non- 
local fluctuations, leads to self-consistent and complete sets of evolution equations 
without back-reaction. This is a completely different outcome in comparison with the 
average from Buchert's formalism ( J 4) p [rf ) ] (notice the subindex p ). Applying Buchert's 
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average to the Raychaudhuri equation \%a\ yields a similar equation as (66), but with 
a non-vanishing back-reaction term [411 I4zl 14"5] 

4-7T 

(H)- p = -{'H)l-—(p) p + 2Q p [r b ], 



3 

with : 



Q P [n] = ((H - (H) p ) 2 ) p - ((H - U q f) p ± 0, 



(68) 



while the application of Buchert's average to the fluid flow equations (8c)-(8d) yields 
similar back-reaction terms (the energy balance equation (86) remains with the same 
form). In particular, the Hamiltonian constraint (17) takes the form 



8tt 



{n) 2 p [n] = —(p) P [n) - (K) p [r b ] - Q p [ n }. 



(69) 



However, if we consider averaged scalars (A) p and fluctuations constructed with this 
average in the framework of Buchert's formalism, we do not obtain a complete self- 
consistent set of evolution equations unless we make extra assumptions on the back- 
reaction terms [4"2"1 133"! ISO] . 

Evidently, since invariant scalars like a ab a ab , E ab E ab and a ab E ab are quadratic 
fluctuations with respect to the q-functions p q , H q (which relate to the q-averages 
(p) q and (H) q ), and not with respect to functions associated with (p) p and (H) p , 
quadratic and higher order fluctuations in Buchert's formalism are not equal to these 
invariants, and thus we cannot obtain in this formalism the simple straightforward 
relations between averages of these invariants and statistical moments of p and H. of 
section 6. However, we can still obtain analytic expressions for computing Buchert's 
average for these invariants. Considering that for any domain the q-average and 
Buchert's average of any scalar are related by 

(AF) P [n] 



(A) q [n] = 



(70) 



together with the identities (52a)-(52c) and (53a)-(53c), we can obtain the following 
exact expressions for Buchert's average of quadratic scalar contractions: 

{H 2 F) P {UF) 2 p Cow p {a ab a ab ,F) 



(o- a bO- ab ) p 



(E a bE ab ) p 



32tt 2 



(P 2 F) P {pTf v 



{<J ab E ab ) p = 8tt 



{pUF) p ( P T) P (HT) P 



Cov p (E ab E ab , F) 

Cow p {a ab E ab , 
{?)p 



(71a) 
(716) 
(71c) 



which lead to the following appealing identities which relate the two averages of these 
scalars and their correlation with the weight factor T: 

Cov p (a ab a ab , T) 

(^)p 

„,,., _ Cov p (E ab E ab , F) 



[O-abO- 



\0~abO- 



(E ab E 



ab\ 



(E nb E 



(<TabE 



ab\ 



(o-abE 



a b \ 



Cov p (a ab E 



(72 a) 
(726) 



V 
ab 



(72c) 



where we have omitted the [r b ] symbol to simplify notation and Cov p is the covariance 
correlation moment (42) with respect to Buchert's average. 
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9. A formalism of exact perturbations on a FLRW background. 



It is evident that (as pointed out in previous work [35l|36]) the system (58a)-((58d| 



has the structure of evolution equations of spherical dust perturbations (the S^ A> ) on 
a FLRW background (defined by the A q ). While the non-local perturbations were 
not considered i n the s e refe rence s, th e same resemblance to dust perturbation holds 



(646| and (65a)-(656| involving (A)q[rb], A q {r) and ( r i r b) 



for systems like (64a 

Following the standard methodology [57J HH1 OHO E31 IM]j a perturbation 
formalism linking a "lumpy" spacetime (LTB model) and a homogeneous "background 
spacetime" (a FLRW dust model) in a domain T>[r] can be defined by comparing local 
LTB variables with other LTB objects that define an FLRW "background state" of 
"zero order" variables by means of suitable maps. However, such maps must also deal 
with gauge issues involving a specific time slicing and coordinates. Since the boundary 
of every spherical comoving FLRW region can be mapped into the boundary of a 



domain T>[r] of an LTB model by the Darmois matching conditions (31o)-(316|, and 



bearing in mind that LTB models and dust FLRW universes are both (i) spherically 
symmetric, (ii) have a geodesic 4-velocity and (iii) their full dynamics reduces to scalar 
modes [30, 32 , a perturbation formalism associated with the q-scalars in domains T>[r] 
can be defined rigorously by means of maps between covariant scalars with all gauge 
issues resolved. 



9.1. The perturbation maps. 

Let Ai be an LTB model and X(T>[r\) the set of all covariant scalars in an arbitrary 
domain T>[r] of Ai. Let Ai be a dust FLRW model and X(Ai) the set of covariant 
scalars of Ai: 

The local perturbation map. For every T>[r] in Ai there exists a model 
M such that A q (t, r) = A(t) holds for every A e X(M) and A e X(M) and 
for all t. The following maps 

$ : X(M) -> X(V[r}), A h> $(i) = A q (r) e X(D[r]), (73) 



<*< ) : X(V[r]) -> X(V[r}), A4^' = 



A-$(A) = A-A q = A- {A), 

$(i) A q (A) q 



(74) 



define for every T>[r] a "background state" associated with an FLRW 
cosmology Ai and local exact perturbations of scalars A obtained by 
comparing them with the LTB scalars produced by the map $. 

An analogous scalar perturbation formalism for the non-local fluctuations can be 
defined along the lines of ( 73 1 and ( 74 ) : 

The non-local perturbation map. Let F(X>[r;,]) be the set of all linear 
functionals in an arbitrary fixed comoving domain Z>[r&] C 3 T[t}. For every 
FLRW covariant scalar A € X(A4) and every T>[ri] the following map 

$ NL : X{M) -> Y(V[r b }), A^(A) q [r b ], (75) 

defines a "background state" associated with a FLRW cosmology Ai 
but consisting of the functionals (13). The definition of the non-local 
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perturbation of the scalar A is then the map 8^1 : X{M) x _X"(Z>[r(,]) — > 
X(V[r ]) x Y(T>[r ]) such that 

0,A)»&-*^^-*&=Vm, (76) 

which provides a pointwise comparison between the real valued functions A 
and their associated functionals (A) q [rb) assigned to the whole domain T>[r^\. 

It is important to emphasize that in the definitions above we distinguish between 
the FLRW "background state" and the FLRW "background spacetime" M. As 
commented in section 4.4, the former is defined by the q-scalars: {p q , H q JC q , fl q } 
or {(p)q[ r b\, (H)q[ r b], (fc)q[ r b], (tyq[ r b]}, which satisfy FLRW dynamics and relate to 
the FLRW scalars {p, H, K,, f2} of the background spacetime M. by continuity under 
the Darmois matching conditions at an arbitrary constant r. These q-scalars are (in 
the perturbation maps) the "zero order" variables, with the "first order variables" 
being the local fluid flow scalars A related to the former by the perturbations or 

S (A) 



Notice that the perturbation formalism defined by ( 73 1 and (74 1 is covariant 
because the background variables and the perturbations (zero and first order variables) 
are coordinate independent LTB objects (see sections 5 and 6). In fact, following the 
Stewart-Walker lemma [23 dE HE G>S] , the perturbations and <5nl^ are also gauge 
invariant (even in the usual sense of 66 ) as they vanish in the FLRW spacetime M. 



associated with the background state through ( 73 1 and ( 74 ) 



9.2. LTB models as exact perturbations. 

As opposed to the conventional approach to perturbations, either the traditional 
one with gauge invariant variables [531 [Ml G>5] or the covariant formalism of Ellis 



et al [23EH1I2S]), the perturbations that emerge from (|73)-(|74|) and (|75|)-(|76|) are 
exact (not approximate) quantities, and thus do not lead to some unknown "near 
FLRW" space-time on the basis of a linearization process (though a rigorous linear 
limit can be defined, see section 12). Instead, the {A q , 5^} or the {(^4) g [r ], 5nl } 
with A = p, W, 57, JC express a known class of spacetimes (generic LTB models) as 
exact spherical perturbations on an abstract FLRW background state defined by LTB 
objects: the A q along all domains X>[r] with varying boundary or the (^4) 9 [rf,] at a 
fixed domain. Evidently these perturbation formalisms are special in the sense that 
they preserve the spherical symmetry and the dust source of the FLRW background. 
However, these formalisms are applicable to any space-time compatible with an LTB 
metric in the comoving frame and having an anisotropic fluid source [36) . and can be 
readily generalized for the non-spherical Szekeres dust models [62] . 

10. Confined perturbations, Swiss cheese models and asymptotic 
perturbations. 



In the definitions ( 73 )— ( 74 1 and ( 75 )— ( 76 ) we assumed domains V[r] or V[rb] with r 



and r\, finite. As we show below, the resulting perturbations can also be considered 
in the asymptotic case when r\> — > oo. 
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Figure 1. Local confined fluctuations. The fluctuations D(A) (thick vertical 
lines) follow by pointwise comparison of the scalars A and A q along r in an 
arbitrary time slice 3 7"[t]. Clump and void profiles are displayed in panels (a)- 
(b) and ( c)-(d), the case of a Swiss cheese model matched to a FLRW region 
at r = rb is displayed in panels (b) and (d). Notice how in these latter cases 
D(A) = at r = rj, thus forcing A q (ri) = A(ri,), which explains the "humps" 
and "bags" in the radial profile of A. 



10.1. Perturbations confined in comoving domains. 

As long as r and rj, are finite, we are considering perturbations that are confined 
in arbitrary bounded comoving domains of an LTB model. Depending on whether 
we perform a smooth match with a FLRW spacetime or not we have the following 
possibilities: 

• There is no matching with FLRW (sec panels (a) and ( c) of figures 1, 2 
and 3). The perturbation describes the dynamics of the domain T>[r] or T>[rb] 



with respect to a FLRW background state defined by ([73| or (75), which is 
a fictitious reference FLRW dust model M. (the background spacetime) whose 
scalars A match (via Darmois matching conditions) the zero order variables A q (r) 
or (A) 9 [r fc ]. Notice that: 

— The fictitious reference FLRW dust model M. is necessarily different for 
different domains. 

— First order quantities like A and the gradients A' and A' q {A = p, H, IC, £1) 
are not continuous at the domain boundaries r or r^. See section 4.4 and 
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Figure 2. Non-local confined fluctuations. The panels display (as in figure 
1) clump and void profiles with and without a Swiss cheese matching to a FLRW 
region at r = rj in an arbitrary time slice 3 T[t}. The non-local fluctuations 
D NL (A) (the arrows) compare local values of A(r) with the average (A) ? [r(,] (thick 
line) assigned to the whole domain X>[r"j,] (shaded region). Notice (as in figure 1) 
that T> N i J (A) only vanishes at r = rj, if there is a Swiss cheese matching (panels 
(b) and (d)). 



panels (a) and ( c) of figures 1, 2 and 3. 
• Swiss cheese holes smoothly matched to a FLRW region (see panels 
(b) and (d) of figures f, 2 and 3). The perturbation describes the dynamics of 
a comoving domain Z>[r&] with respect to a background state that in this case 
corresponds to the actual (non-fictitious) FLRW background spacetime matched 



at r = Tb (under conditions (31a)-(31c|). Notice that 



The resulting spacetime is a compound Swiss cheese configuration consisting 
of an LTB section (confined in 2?[rj,]) described by zero order variables A q 
and their perturbations 5^ or and the FLRW background spacetime 
that extends for r > r&. The FLRW dust model M. is necessarily different 
for different domains 2?[rh]. 

Darmois matching conditions only require the zero order quantities A q and 
(j4) g [r&] to be continuous at r&, though it is always possible to demand (as 
extra conditions) that the local scalars A and the gradients A' and A' q (which 
are first order quantities related to 5^ A > via (35)) be also continuous at 
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rb. These extra conditions imply that the fluctuations and perturbations 
themselves vanish at r b (see section 4.4 and panels b and d of figures 1, 2 
and 3). 

10.2. Asymptotic perturbations. 



The application of the perturbation maps ( 73 )-( 74 ) and ( 75 )-( 76 ) to the asymptotic 



limit r — > oo or r b — > oo depends on the convergence of LTB scalars and their 
associated q-scalars along radial rays (i.e. curves x a = 15%, where t = J g rr <ir is 
the proper radial length), which are spacelike geodesies of the metric ([lj and of the 
slices 3 T[t] with metric h ab . As shown in |39) . t(r) is a positive monotonic function in 
regular LTB models, and thus the asymptotic limit r — > oo corresponds to £ — > oo. 

The radial asymptotic behavior of LTB models is determined by the following 
asymptotic limits for the scalars A = p, %, fC, 0: P~| 

( A x (t) ^ 0, Asymptotically FLRW 
lim A(r) = lim A q (r) = lim (A) q [r b ] = I (77) 

[ 0, Asymptotically Minkowski. 

These limits correspond to q-averages for increasingly large domains up to the 
situation in which T>[r] or T>[ri,} become the whole slice 3 T[t], and thus, we can 
speak of the q-average of a whole LTB model (instead of the q-average of confined 
domains of an LTB model). We look at the perturbations for asymptotically FLRW 
and Minkowski models separately below. 



Asymptotically FLRW models (sec figure 4). As a consequence of (74) and 



(77) local perturbations vanish asymptotically in these models: 

lim S^(r) =0. (78) 



However, (76) and &771 lead for non-local perturbations to the following non- 



trivial nonzero asymptotic limit: 

AA), x Air) -A (t) Air) -Ait) „ . 

lim 6£\r,r b ) = W °° V ' = V ' W , r finite, 79 
r »^°° A^it) A(t) 

which keeps their interpretation of a "contrast" between local values A(r) and the 
domain average, with the difference that the average now corresponds to a domain 
that encompasses the whole slice 3 T[t] and coincides with the corresponding 
scalars A^ (t) = A(t) of an asymptotic FLRW background state defined by 



(75 1 . As a consequence, we can rigorously state that the q-average of every LTB 
model of this class is exactly the FLRW background spacetime Ai. Asymptotic 
fluctuations and perturbations for a model converging to FLRW are illustrated 
by figure 4. 



Asymptotically Minkowski models. By looking at (76) it is evident that 
non-local perturbations cannot be defined in the asymptotic limit for these 
LTB models, since in general A(r) ^ for finite r and thus the ratio A/(A) q [rb] 
diverges as (A) q [ri 1 } — > in the limit r b — > oo. However, local perturbations 5^ A ' 
defined by (l74|) are not affected by the radial convergence to a Minkowski vacuum 



+ LTB models can only converge in the asymptotic radial regime to FLRW dust models with zero or 
negative spatial curvature, not positive. Models convergent to Minkowski can converge to a section 
of Minkowski parametrized by Milne coordinates or by non-standard curvilinear coordinates. In the 
Milne case p q , f2 9 — > but H q , K q tend to nonzero values (see 1391 ). 



Averaging and exact perturbations in LTB dust models 



24 



because, as shown in [39], both A — > and A q — > hold as r — > oo but their ratio 
A/A q is finite in this limit, and thus the 8^ A > in (74 1 are well defined. The only 



difference with asymptotically FLRW models is that the 5^ A ' tend (in general) to 
nonzero constant values as r — ► oo (see [39 ). 

11. Comparison between local and non local perturbations. 

Since the perturbations 5^ and the #nl are different objects, they provide a different 
measure of the inhomogeneity or "deviation" of an LTB model from a FLRW Universe: 

• The 5^ measure this deviation through local comparison of A(r) with the 
characteristic FLRW scalar A q (r) = (A) q [r] at every r, and is related (via (35 1) 



to radial gradients of A q and A and (via (49a)-(496|) to the ratio of Weyl to Ricci 



curvature and anisotropic to isotropic expansion. See figures 1 and 3. 

• The (5n L ^ measure the deviation from FLRW though the notion of the "contrast" 
between local values of A(r) inside an arbitrary but fixed comoving domain 
(0 < r < r b ) and a value (A) g [rf ) ] that can be associated to the equivalent FLRW 
scalar A(t) = (A) g [rb\ that characterizes the whole domain. In the asymptotic 
limit r& — > oo so that A(t) — A^ (t), these perturbations provide the contrast with 
respect to the FLRW spacetime to which the LTB model converges asymptotically. 
See figures 2 and 3. 

We examine below how these differences relate to the shape of the radial profile of 
a scalar at an arbitrary 3 T[t] and the signs and magnitudes ("amplitude") of the 
perturbations. 

11.1. Signs of the perturbations vs type of radial profiles. 

Consider a scalar A whose radial profile is monotonic throughout a given domain [^] 



Then, following (23 1, the clump/void profiles arc defined by 

Clump profile A' < =*> A' q < A < A qi (80a) 

Void profile A' > => A' q > A > A q , (806) 

where A(Q) = A q (0) — A c and A'(0) — A'(0) — and we have assumed absence of 
shell crossings and singular layers (hence R > holds everywhere or it has a common 
same order zero with A' and A' [3S1 SO]). The relation between the signs of the 



perturbations and the profile type follows from (351 (it can also be seen in figures 1, 
2 and 3): 

• Local perturbations: 

clump profile : 5^ < 0, 
sign(^) = sig n(4) ^ { (81) 

void profile : S (A) > 0. 

* Such domains always exist for n finite, even if A is not monotonic for values r > r^. The nature 
and evolution of radial profiles of LTB scalars were examined extensively in reference 1401 . In general, 
density void profiles are fully compatible with regularity conditions (absence of shell crossings) for 
hyperbolic models, but not for elliptic or parabolic ones. Also, for models with a non-simultaneous 
big bang (t' bb ^ 0, nonzero decaying modes) density void profiles always emerge from initial clump 
profiles after a transition ("profile inversion"), whereas a density void profile exists for all the time 
evolution only in models with a simultaneous bang (zero decaying modes). A profile inversion of 
the expansion scalar H necessarily occurs in elliptic models but not in hyperbolic ones (see |40| for 
further detail). 
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• Non local perturbations: 

Finite domains T>[r b ] and Swiss cheese holes 



The sign relations in (80 a I and 
(80&I imply that there always exists a value r = y < r b such that A(y) = A q (r b ) = 
(A) q [r b ], and thus we have from (39): 5i^\y,r b ) = 0. Since A c > (A) q [r b ] holds 



in clump profiles and A c < (A) q [r b ] in void profiles, then dm? changes sign as 
follows (see panels (b) and (d) of figures 2 and 3): 



Clump: 4l ) > for < f < y, 6£' < for y < f < r b , (82a) 
Void: 4l 3 < for < f < y, 5$ > for y < f < r b , (826) 

Asymptotic background (see figure 4). In the limit r b —¥ oo we have A — > 

3 T[i], hence for all r: 

(83 a) 
(836) 
836| clearly 



A Ml (A) q [r b ] -> and V[r b ] 
Clump: A > A^ 
Void: A < A„ 



> 0, 



sii ] < o, 



Comparing the sign relations in (81 ) with those in (82a)-(826) and (83a 
indicates that non-local perturbations S^t have the expected intuitive signs that we 
tend to infer for perturbations through the familiar notion of a "density contrast" with 
respect to a fiducial (or background) FLRW density value p: an over-density (clump) 
is a positive perturbation because local density is larger than p, while an under-density 
(void) is a negative perturbation because p is smaller than local density, which we can 
now identify with (p) q [r b ] or its asymptotic limit In fact, as shown in figures 2, 3 
and 4, this intuitive notion of a "contrast" can be applied to all covariant scalars, not 
just the density. On the other hand, the "gradient" based local perturbations 8^ 
have the opposite sign to the intuitively expected one and thus are more abstract and 
counter-intuitive (see figures 3 and 4). 



11.2. Amplitudes. 



The difference between the perturbations 5^ and S Nl , can also be understood in 
terms of the notion of "amplitude" often introduced in the context of simple intuitive 
Newtonian [ST] and relativistic [31131 [TO] perturbations. In order to illustrate this point 
we consider an asymptotic perturbation in an LTB model converging radially to FLRW 
(as in figure 4), so that for all A we have along the slices 3 T[t] the asymptotically 
FLRW behavior of @: 

lim A(r b ) = lim (A) q [r b ] = A^t) = A(t) £ 0, (84) 

where A(t) is the A scalar for the asymptotic FLRW state (see figure 4). Assuming 
for simplicity a monotonic profile, then without loosing generality we can express any 
A in any complete slice 3 T[t] (not intersecting a singularity) in terms of a "contrast 
amplitude" 4>A as 



A = A^ [1 + (f> A f(r)] , with A =4f(O,oo) 



lim 



A r 



{A) q [r b ] 



A 



(85) 

where A c = A(t, 0) = A q (t, 0) and / is a smooth non-negative function satisfying 
fc = /(0) = 1 and /'(0) = 0, as well as /(r), f'(r) — > as r — > oo. Since we can 
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Figure 3. Confined perturbations. The panels display the local and non- 
local perturbations (<5^^ and <5nl ) f° r the radial profiles displayed in figures 1 
and 2. Notice that both types of perturbations vanish at r = when there is a 
Swiss cheese matching (panels (b) and (d)). Notice also that 5^ A ^ and have 
opposite signs around the center. 



always choose at any 3 T[t] the radial coordinate as R = r, so that R 2 R'dr = r 2 dr, 
the form (85) implies from (13) and (18) 

A q (r)=A oo [l + <j> A f q (r)}, f q (r) = {f) q [r] 



J r f 2 df ' 



(86) 



while (|34|), (|39|) and (|63|) yield 

/ 



(/ - U 



(87) 



l + 4>Afq ' 

which illustrate the following basic features: 

• The contrast amplitude 4> A obeys a simple linear relation with the more intuitive 
non-local perturbations but its relation with local perturbations 5^ is 
non-linear. This is consistent with the fact that the latter perturbations are less 
intuitive than the former. 

• Both perturbations vanish in the domain boundary r — > oo where the 
correspondence with the FLRW background occurs, but their behavior at the 
center is different: 5^t' (0, oo) = 4> A ^ (in general) while 6^(0) = (see figure 
4). 
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Figure 4. Asymptotic fluctuations and perturbations. Panels (a) and (b) 
correspond to a clump profile, while ( c) and (d) depict a void profile. Since the 
domain encompasses a whole time slice 3 T[t], non-local fluctuations compare local 
values of A(r) with the global average (A) 9 [r;,] — > A^. in the limit — ¥ oo, which 
is equal to the corresponding scalar A(t) of the FLRW asymptotic state. Non— 
local perturbations <5^ L yield the "contrast" between local values of A and this 
average. The amplitude of the perturbation is (ft a, while ty A depicts the maximal 
value of the local perturbation jC 1 *). Notice that S^j} is positive/negative for the 
clump/void, while 5^ A ' has opposite signs. Also, both perturbations vanish as 
r — > 00. 



• The relation between the signs of the perturbations and the type of profile can 
also be given in terms of the amplitude: 

Clump profile: A c > A x , 4> A >0, 5$ > but 6 {A) < 0, (88a) 

Void profile: A c <A ao , (j) A <0, 5$ < but 5 {A) > 0. (886) 

wher e we u sed the fact that (by construction) /' < so that the sign relations 
(80a|-(806) imply / — f q < and the sign of 6^ A ' is the opposite of the sign of 



Evidently, the convention whereby a clump or void (i.e. overdensity or underdensity 
when A = p) respectively correspond to positive or negative contrast follows from the 
sign of the amplitude contrast embodied in <5nl\ and is more intuitive and easier to 
follow than the complicated relation between radial gradients A' q associated with S 1 -^. 

However, in spite of their differences in signs when referred to clump/ void 
radial profiles, the deviation form homogeneity (FLRW conditions) can be traced 
effectively with both types of perturbations, S^} and 5^ , because there is a consistent 
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Figure 5. Magnitudes of local and non-local perturbations. The graph 
depicts the monotonic non-linear relation between the amplitude, </>a , of non-local 
perturbations and the maximal value of the local perturbation, tpA (see figure 
4), which shows that both perturbations provide a self-consistent qualitatively 
analogous description of the deviation from FLRW. Only for small amplitudes 
\4'a\ *C 1 the relation is linear. 



monotonic relation between their magnitudes (their maximal/minimal value in a given 
domain) and the amplitude <f>A- For the non-local perturbations 8^} the amplitude 
4> A coincides with their extremal value, which occurs at the center r = (for domains 
with monotonic /. For local perturbations the relation between (j)A and the extrema 
of 8^ (denoted by ip A , see figure 4) is more complicated, as the latter may occur for 
different values of r for different slices 3 T[t] and depends on the choice of /. Since 
a result for a general / is hard to find, we consider the special form / = (1 + r™) -1 
where n > (the results are qualitatively analogous for all other forms). The result, 
shown in figure 5, illustrates the non-linear monotonic relation between ip A and 4>a- 
Hence, both types of perturbations provide a consistent estimate of the deviation 
from homogeneity of the models in which larger values of \4>a\ correspond to larger 
magnitudes of 8^ and 8^. 



12. Linear limit. 



While the <5 (j4) and the ^ni? are exact quantities that need not be "small' and do 
not comply with linear evolution equations, we can expect (intuitively) that they 
should somehow reduce to linear dust perturbations (in the comoving gauge) when 
their magnitudes (i.e. amplitudes) are small. A more rigorous way to examine their 
connection to linear perturbations of dust sources follows by constructing second order 
equations for 8^ and 8^: 



• Local per turb a tion s. We differentiate both sides of (58 a I and use the remaining 
equations (586)-(58d| to eliminate all derivatives except 8^ and S^ p \ leading to: 



2 [j^] 2 



1 + 8(p) 

• Non local pert urbations. 



2U q S {p) -^ Pq 8 (p) (l + jW) 



0. 



(89) 



We follow the same procedure as above: 



differentiate (|64a|) and use (|646|-(|65fe|) to rewrite it in factors of 8^} and Sst ■ 
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Considering asymptotic backgrounds (rf, 
2 [8 



■'NL 



1 



dp) 



4(6H o 



m q )s 



(p) 

NL 



oo ) yields: 



47rp 0O 4 P L ) + 6(K 9 -K DO ) i 



1 + 8 



= 0, 
(90) 



which holds for confined domains by replacing p x , with (p)q[rb\, (^)g[ r f>]- 

We remark that both (89) and (90) are exact non-linear equations for 5^ or 8^ 
(equations similar to (89) have been obtained in [551 for Szekeres models). Also, 
notice that ( 89 ) and ( 90 1 coincide as r — > oo for asymptotic perturbations and as 
r — > r^ for confined domains T)[r^. 

In general there is no reason to assume that the amplitudes §a or their time 
derivatives 4>a are small in a generic LTB model subjected to non-linear evolution 
equations, but if in a given range of t we have \4>a\ ^ 1 then the relation between 
both 8^ and 8^} with <pA should bec ome similar. This can be seen by expanding for 
\4>a\ <^ 1 hi the non-linear relation in (87) between the 8^ and the amplitude 4>a'- 

<5 (A) «(/-/,) 0a + O(0i), (91) 
so that the 8^ become linear on cf>A (as the 8^/), and thus the magnitudes of 8^ 
and 8^ are both proportional to <Pa because < / < 1 and the fluctuations f — f q are 
bounded. As a consequence, if \4>a\ ^ 1 both perturbations are of the same order in 



4> A , and thus both comply with \8^\ -C 1 and \8^} 
of 5^ must be the same as that of 8^} 



<C 1 and the time evolution 
at leading order in <p A . In particular, 
if we take A = p, H, the linearity conditions \<j> | -C 1 and \4> H \ <C 1 imply that 

\8^\ « 1, \Sf$\ « 1 and \8W\ « 1, |4? 3 | « 1 hold, and thus \8^\ « 1, |^| « 1 

64a[ )). Therefore, we also have p q « p^ « p and 
% (the same relations hold for finite domains by replacing p^ , H x 



also holds (from (58c 



with (p) q [ro], (H) q [r ]). Considering all these implications, if \4>a\ -C 1 the second 
order evolution equations ( 89 1 and ( 90 ) becomes at order <fi A 



2 ^oo 5W ~ 4np x = 0, (92) 
which is formally identical to the evolution equation for gauge invariant linear density 
perturbations of a dust source around a FLRW background characterized by p^ , 
(or (p)q[rb\, (H)q[rb] for bounded domains) in the comoving gauge [67l[70], which for 
dust is a synchronous gauge as well. 



13. Gravitational entropy. 



A gravitational entropy functional has been proposed by Buchert et al. [42, 57, 581159). 
as an application of the Kullback-Lciblcr functional of Information Theory to 
spacetimes with inhomogeneous dust sources. While Buchert and coworkers defined 
this functional with the standard proper volume average (i.e. (13) with T = 1), as 



this is the scalar average of Buchert 's formalism, we consider here the same functional 
with the weighed q-average (13) with T ^ 1. Hence, our approach will yield different 
theoretical connections. 

Following Buchert et al, we can define the Kullback-Leibler entropy functional in 
terms of the q-average over a domain T>[ro] as 



S - 5 cq = 70 / p In 
>T>[r b ] 



{p)q[n\ 



FAV P = 7o ( p In 



(p)q[n] 



[n]V g (r b ), 



(93) 
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with: V q = [ FdV p = ^-R\r b ), (94) 

JV[r b } -3 

where 5 cq = S ccl (r b ), and we have introduced the constant 70 = k B /(mc 2 ), with k B 
being the Boltzmann's and m a particle mass, so that S — S B€l has units of entropy. 
If we define the entropy current associated with (93) as S a — p(S — S ccl )u a , then 
by considering the rest mass density conservation law (86 1 we obtain the following 
entropy balance law: 

V a S a > =► S>0, (95) 

which allows us to identify the inhomogeneity measured by p(r) ^ (p)q[ r b\ &s the 
source of gravitational entropy production (S > 0) in arbitrary domains P[r&] of LTB 
models, with equilibrium states complying with S = 0, S = 5 cq and corresponding to: 

• FLRW models as the homogeneous subset (particular case) of LTB models for 
which p(r) = (p) q [r b ] holds for all r, r b , t. 

• specific boundaries of generic LTB models, such as: 

— The center worldline (in models admitting a regular center) corresponding 
to the case r b = 0, since p(t, 0) = p q (t, 0) holds for all t [39l H0j (56] . 

— Asymptotic boundary r — > 00 along radial rays of time slices 3 T[i] of LTB 
models radially converging to a FLRW state (see section 10.2 and [33]), for 
which both p and p q tend to the same value p — p(t) in this limit (see figure 
4). 

— The boundary r = r b (for all t) of Swiss cheese holes in domains Z>[r&] 
smoothly matched with a FLRW region when the extra conditions ( 32 ) are 
imposed, so that p(rt) = p q (rb) holds (see sections 4.4 and 10.1, and panels 
(b) and (d) of figures 1, 2 and 3). Evidently S = holds for all t in the 
FLRW region r > r b . 

We remark that in all the equilibrium states listed above we have £W = 4l = 0, and 
thus entropy production is closely related with nonzero density perturbations. For the 
remaining of this section we will assume domains such that 5^0 holds. 

In order to evaluate S, we apply the time derivative commutation rule ( 22 ) to 
( 93 1 , leading after some algebraic manipulation to the same relation between S and 
the non-commutativity of the time derivative and the average found in |57l 1581 159j : 

' S ' ( ' 7 ' ! ={p) q [n]-{pUnl (96) 



ioV q (n) 



which, with the help of (22) and (|42j), can be related to the negative statistical 
correlation of fluctuations of p and "H, given by the covariance momentum with respect 
to the q-averages of these variables in a domain 2?[rt]: 

~ = -3 7o [{ P H) q - ( P ) q (H) q ] = ~3 l0 Cov q (p,H) = -37o(D(p)D(W)>, > 0, (97) 

Vq 

so that: 

Cov q (p,n) = (D(H)-D(p)) q [r b ] < => S(r b ) > 0, (98) 

where we have applied ( [42] ) to replace non-local fluctuations with local ones (see proof 
in Appendix C). Considering ([5])-([6]), (44), (46c) and (52c), condition (98) can be given 
in terms of the q-average of a scalar invariant by: 

{(TabE ab ) q [r b ] <0 S(r b ) >0, (99) 
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which is a very elegant way to connect ( 93 ) with an unequivocal and completely 



coordinate independent marker ol inhomogeneity, as it contains contributions from 
density and velocity fluctuations 

in equilibrium states implies 8^ p ' 



condition S 



Notice that, as consequence of (97)-(|99|) 

km 



?0) 



the 



0, and 



thus perturbations of all q-scalars necessarily vanish in these states. 

As pointed out by Buchert et al in. [571 EHl EH] (who also obtained the relations (97) 
and (98) but not (99)), a negative correlation of the density and Hubble fluctuations 
seems to be consistent with the intuitive behavior of gravitational clustering processes 
in structure formation scenarios: over/under dense regions tend to contract/expand 
as density increases/decreases. However, this intuitive behavior occurs in the time 
evolution, not (necessarily) in domains along time slices where averages are computed. 
Therefore, as much as ( |98[ ) is an intuitively appealing and elegant proposition (in view 
of (99)), it is not a priori evident that we will obtain the desired negative sign of 



D( j o)D('H) when actually computing the fluctuations. Since verifying (98 1 in arbitrary 



domains of generic inhomogeneous LTB models is a comprehensive task that is beyond 
the scope of this paper, and thus will be accomplished in a future separate article [71] , 
we provide in this section a guideline on how this verification can be accomplished, as 
well as some preliminary results (listed in the following section). 

Testing ( 98 ) involves computing radial integrals of quantities whose explicit radial 



dependence is not known (and thus must be done numerically). However, we can infer 
the fulfillment of ( 98 1 by means of the following sufficient (but not necessary) condition 
that can be evaluated at the boundary of every domain r = r a : 

D(ft)D(p) = [H(r b ) - U q {r b )][p{r b ) - p q {r b ) ] < => S(r b ) > 0, (100) 



where we used ( 42 ) (see its proof in ( C.6 )— ( C.7 ) ) to be able to consider the fluctuations 



evaluated at the boundary r = r b of a fixed but arbitrary domain. Considering ( 23 ) 



and ( 24 ) , the product of these fluctuations is 
p' q (r b )H' q (r b ) 



DCW)D(p) 



[V'(r 6 )/V'(r 6 )P 



/ P 'v q df / n'Vgdr, 

Jo Jo 



(101) 



and thus condition ( 100 1 can be given directly in terms of the correlation of the radial 



gradients of the density and the Hubble expansion scalar: |] 

p' q {n)H' q {n) <o => s(r b )>o, (102) 

or (from (35)) in terms of the local perturbations S^ p \ 5^'; 

S ip) (r b )S in) (r b ) < =>> S{r b )>0. (103) 

Either one of ( 102 ) or ( |103[ ) allows us to test (at least qualitatively) the fulfillment 
of S > in arbitrary domains of generic LTB models by means of the analytic and 
qualitative results derived in [SHI [3U SO]- We summarize in the next section some 
preliminary results (further results and more detail will be provided in |71j). 



Besides the balance law (95), it is useful to compute ti me c oncavit y /convexity of 



S given by the sign of S, which follows readily by applying (22) to (97) 
S 



97o 



Va 



p[B(H)] 



; H q T>( P )T>(H) 



4ir 



D(p)]< 



t) The condition p'(r)W(r) < is a lso a sufficient condition if it holds for t he ful l inte gration range 
< r < r - ;, in the integrals in | |101| . On the other hand, conditions ( |102[ l and | |103| only need to 
be evaluated at tv Notice that for monotonic profiles p'(r b ) < and p'(rt) > respectively imply 
Pq( r b) < and p'„(r b ) > 0, but the converses are false, see |40| . 
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Figure 6. Time slices and the expanding and collapsing singularities. 

The figure displays time slices ti, t2, t^, (long dashed horizontal lines) in the 
(t, r) plane and the worldline of a comoving layer r = constant (short dashed 
vertical line) evolving from the a non-simultaneous big bang (thick lower curve 
tbb{r)) towards the collapsing singularity (thick upper curve t co n('>')) in an elliptic 
model. For t\ , f 2 we have necessarily H' q < and p ' n < 0, whereas we have H' q > 
and p' < for ti, t2- Therefore, conditions ( 102 I and 1 103 1 hold for slices near 



t co ii but fails for slices near t bb . The latter failure also occurs in hyperbolic models 
with non-simultaneous big bang that lack a collapsing singularity. 



370 
4 



P o- ab a ab + ^-H q a ab E ab + -^-E ab E ab 



1 

in' 



(104) 



where we used (8a)-(86|, (44|, ( 576 ) — ( 57fo I and (58a)-(586) to eliminate the derivatives 



p, p q , H, H q and S in terms of fluctuations D(p), D("H). 



14. Conclusion and summary. 



We have considered the q-average (the scalar average defined by (13)) in spherical 



comoving domains T>[r] for the study of generic LTB models and their properties. 
This average generalizes the standard proper volume average of Buchert's formalism 
by the presence of a non-trivial weight factor, J 7 , that is suitable for these models 
(Buchert's average corresponds to the case T = 1). By applying the q-average to 
LTB covariant scalars A, we obtained new dynamical variables (q-scalars) that can 
be defined as functionals (A)q[r b ] (the average is assigned to a whole fixed but arbitrary 
domain) or as functions A q (q-functions that depend on the average for domains with 
varying boundaries). The q-scalars are covariant, as they are expressible in terms of 
curvature invariants (see section 6), and also have very special and useful properties. 
Specifically, the q-scalars that are common with FLRW models (A q and (A) q [r b \ for 
A = p, H, IC, ft) identically satisfy FLRW fluid flow dynamics and their analytic forms 
mimic FLRW expressions that are commonly used in LTB void models that probe the 
possibility of explaining observations without resorting to dark energy (section 4.3). 
Also, the q-scalars determine (via Darmois matching conditions) a unique FLRW 
background state (section 4.4) that can be either a fictitious reference FLRW model 
that is different for each domain, or an actual background FLRW spacetime if a domain 
T)[r b ] is smoothly matched to a FLRW region. Fluctuations and perturbations of the 
q-scalars can be defined by comparing them with non-averaged (i.e. local) covariant 
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scalars. These fluctuations and perturbations are local (D(A), S^ A ') or non-local 
(D nl (j4), <5nl'), depending on whether they are, respectively, defined for q-scalars 
that are functions or functionals. The following results are worth highlighting: 

• All proper tensors (curvature and kinematic) characteristic of LTB models are 
expressible by linear combinations of g a b, u b and the common eigenframe of the 
shear and electric Weyl tensors (e£), with coefficients given by the q-scalars 
A — p q , Hg, K, q , fig and their relative local fluctuations (local perturbations) 
§(p) ^ S^ n \ 5^, 6( n \ In particular, all curvature tensors are expressible in terms 
of these linear combinations with coefficients given by the Ricci scalar 1Z and its 
local fluctuations D(7\L) (section 6). 

• Quadratic local fluctuations of p and % are equal to invariants that follow from 
contractions of the tensors involved in the dynamics of LTB models: the Ricci, 
electric Weyl and shear tensors. The q-averages of these fluctuations (local and 
non-local) are equal to statistical moments (variance, correlation and standard 
deviation) of p and T-L (see section 6). 

• Local perturbations 8^ ' provide an elegant and invariant characterization of 
inhomogeneity in terms of the ratio of Weyl to Ricci curvature: 4 , 2 /^i where 
$2 is the on ly nonzero Newman-Penrose conformal invariant and 1Z is the Ricci 
scalar, whereas relative Hubble scalar fluctuations 5^ do so in terms of the 
ratio of anisotropic vs. isotropic expansion: where £ is the eigenvalue of 
the shear tensor. All q-scalars can be given in terms of these scalar invariants 
(section 6). 

• The evolution equations for the q-scalars and their perturbations (local and non- 
local) completely determine the dynamics of the models, and thus provide an 
alternative to the use of analytic solutions in the study and applications of the 
models. Although the resulting systems involve PDE's, they contain only time 
derivatives and the constraints are purely algebraic in nature and preserved in 
time. Hence, these equations can be treated effectively as ODE's (section 7). 

• The basic difference with the standard proper volume average of Buchert's 
formalism is the absence of back-reaction correlation terms in the evolution 
equations associated with q-averages and their perturbations. Also, the relation 
between the averages of quadratic and higher order fluctuations and averages of 
scalar invariants and statistical moments of p and % does not occur with Buchert's 
average. We derive exact expressions for the Buchert's average of quadratic scalar 
invariants (see section 8). 

• Since the q-variables behave effectively as FLRW scalars and their perturbations 
convey for every domain the (exact) deviation from FLRW dynamics, we can 
rigorously re-interpret LTB dynamics as the dynamics of exact spherical dust 
perturbations on a FLRW background defined by the q-scalars. As mentioned 
before, this background state is an abstract reference background defined at every 
domain by the continuity of the q-scalars under Darmois matching conditions 



(31a|-(31c). It is an actual FLRW spacetime only if a smooth matching with a 
FLRW region is considered in the context of "Swiss cheese" models (section 9). 

Both local and non-local fluctuations and perturbations can be either confined in 
a given domain T>[rb\ (with or without assuming a "Swiss cheese" configuration 
through a matching with a FLRW region, see figures 1, 2 and 3), or an asymptotic 
perturbation for the case when V[rb] becomes the whole time slice 3 T[t] in the 
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limit rb — > oo (see figure 4). The choice between confined and asymptotic 
perturbations depends on the boundary conditions of the specific problem or 
application that we may work out with LTB models. Confined perturbations 
(local or non-local) without a matching with FLRW are practical because they 
are easily applicable to study the inhomogeneity pattern of any generic LTB 
configuration, but are more abstract because the FLRW background state is 
a fictitious reference spacetime that changes for each domain. Swiss cheese 
configurations are more intuitive because the background state is a single (and 
non-fictitious) FLRW spacetime. They are useful if we wish to describe various 
inhomogeneous regions, as in Swiss cheese void models used to fit observations, 
but as as illustrated by panels (b) and (d) of figures 1, 2 and 3, their intuitiveness is 
offset by the artificial need to impose continuity of the perturbations (which must 
vanish in a non-ficticiuos background at finite radius) . Asymptotic perturbations 
are more natural, and also practical if the radial profile of the scalars rapidly 
converges to a given set of FLRW values. 

• The local perturbations provide a measure of inhomogeneity through the 
local magnitude of the radial gradients of covariant scalars, while non-local 
perturbations do so through the familiar notion of the "contrast" between local 
values of these scalars and a reference FLRW value defined as a q-average for 
each fixed whole domain, or for the whole slice in the asymptotic limit. Even 
if their evolution equations are more complicated, non-local perturbations are 
more intuitive than local ones, as they allow us to associate over and under 
densities with (respectively) positive or negative density amplitudes, while local 
perturbations yield the opposite sign. We have illustrated these points through 
the notion of the "amplitude" of the perturbations (see sections 10 and 11 and 
figure 5). 

• The definition of non-local perturbations 8^} by the maps ( 75 1 and ( 76 1 provide 
(through the q-average and its relation to kinematic and curvature invariants) 
a rigorous and coordinate independent interpretation for simple examples of 
perturbations that are based on the notion of a contrast with respect to a FLRW 
background. These simple contrast perturbations are frequently introduced as 
ansatzes in many text-book and articles, for example, in the context of simple 
Newtonian models of structure formation (the "top hat" or "spherical collapse 
model" [67] ) but also in linear perturbations |70j and in astrophysical applications 
of LTB models (see examples in [3| l4| [10] ) . 

• Both local and non-local perturbations yield in the linear limit the familiar dust 
perturbations of linear theory in the comoving (or isochronous) gauge (section 
12). 

As a quick application, we examined in section 13 the definition of gravitational 
entropy based on the Liebler-Kullback functional from Information Theory, which has 
been proposed by Buchert and coworkers [HI \57\ [551 [FJ5] • Considering this functional 
as given in terms of q-averages yields the same conjecture obtained in these references, 
namely: a relation should exist between positive entropy production (S > 0) and a 
negative correlation between fluctuations of p and T~L. However, testing this conjecture 
is easier if we use q-averages, since: S is proportional to the q-average of the scalar 
invariant a a bE ab , and also (p) q and (H) q are expressible in closed analytic forms and 
the evolution equations for their fluctuations and perturbations (sections 7 and 8) 
do not involve complicated back-reaction terms. Hence, we were able to obtain the 
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sufficient conditions for S > in terms of closed analytic expressions, such as ( 102 1 



and ( 103 1, that can be tested in arbitrary domains of generic LTB models (a task that 
will be accomplished in future work [7T]). Meanwhile, as "food for thought" we list 
the following preliminary results: 



Entropy production is positive. Conditions ( 102 ) and ( 103 1 hold in the following 

cases: 



Late time evolution of elliptic models (see (A.3)) in time slices near the 
collapsing singularity. This result can be justified by qualitative arguments 
(see figure 6). 



Asymptotic time evolution (t —> oo for all r) of hyperbolic models (see (A. 2)). 
This result follows from the fact that 5^ — > 8^/2 holds in this limit 
for these models (because Q, q — > 0), and thus from table 1 of [ID] we have 
SkPifiW < for all initial conditions complying with regularity. 



Entropy production is negative. Conditions ( 102 1 and ( 103 1 arc violated in the 
following cases: 



• Parabolic models (spatially flat models K. = 0, Q q = 1, see (A.l)). Since 
§(H) = J(p)/2 holds for all t, and regularity conditions require — 1 < S^ 1 ' < 
(see [3SII0]), then 5^5^ > holds for all t. 

• Early time evolution of elliptic and hyperbolic models with a non- 
simultaneous big bang. This result can be justified by qualitative arguments 
(see figure 6). 

These results clearly indicate an important connection between the failure to fulfill a 
positive entropy production and specific conditions in which the density decaying 
modes dominate over the growing modes (parabolic models and near a non- 
simultaneous big bang [5J |31 HI EUJ 133)- This is an important theoretical feature 
worth a comprehensive examination (further detail is provided in [71]). 

As a comparison with our approach and preliminary results, a recent perturbative 
study [59] of the Liebler-Kullback functional in the context of Buchert's average 
yields S only in terms of the average (Buchert's average) of the invariant scalar 
C a bcd,C ahcd (which is 8E ab E ab in LTB models), hence it excludes the contribution 
to inhomogeneity from the shear tensor (a spatial gradient of the velocity in the 
conformal Newtonian gauge) that comes from the scalar invariant a a t,E ah . Buchert 
and coworkers have not yet applied the Leibler-Kullback functional to generic LTB 
models in a non-perturbative manner (only a very simple example of an LTB model 
was considered in [55] , but this model exhibits shell crossings) . Although they obtained 



conditions that are analogous to (98), their main interest has been to provide a 
theoretical connection between the growth of structure complexity associated with 
S > and the dynamical implications of the back-reaction terms that appear in 
Buchert's formalism. Since the back-reaction terms of Buchert's formalism vanish 
for q-averages, our results prov ide a different theor etical perspective, namely: the 



connection between S > in (97) and S in (104) associated with the growth of 
structure complexity and the statistical moments that arise from the q-average of 
quadratic invariant scalars (E a i,E ab , a a bE ab , o- a \ 3 a ah ) that vanish for FLRW models. 

It is important to remark that the notion of a gravitational entropy connected to 
inhomogeneity (marked by invariant scalars of the Wcyl tensor) bears a close relation 
to the concept of the "arrow of time" originally suggested by Penrose [72] , and further 
explored by a number of authors [731 IZH [7SJ, [73 |77] . The Weyl tensor scalar that 
emerges in this literature is C a bcdC abcd (proportional to E ab E ab in LTB models) and 
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no attempt is made in these papers to connect it with averages and fluctuations. 
As a contrast, the Leibler-Kullback functional is based on an averaging formalism 
that yields a direct connection between the growth of structure (inhomogeneity) and 
its associated fluctuations, which in the case of the q-average, are directly related 
to invariant contractions of a a b and E ab , as well as with fluctuations and statistical 
variance and covariance moments of both p and % (which is a fluid velocity) . Evidently, 
we have only considered this functional with the q-average applied on LTB models, 
and thus, we cannot exclude the possibility that contradictory results may happen 
with Buchert's average or in other types of models (as for example with the "arrow of 
time" , see 75 ) , but the work we have developed here may provide a guideline on how 
to explore the notion of the arrow of time and its relation with the Liebler-Kullback 
functional in more general models. 

While the q-average formalism cannot be used to study and understand the 
dynamical implications of back-reaction, it has resulted very helpful for understanding 
the properties and evolution of LTB models [3^1 EH], including the computation in 
these models of back-reaction terms and verification of the existence of "effective" 
acceleration in the context of Buchert's formalism [56]. Moreover, it does provide 
through its fluctuations and perturbations interesting connections between averaging, 
perturbation theory, invariant scalars and statistical correlations of p and "H, which 
signals a valuable theoretical insight on how the averaging process should work 
in any generic solution of Einstein's equations (at least in LRS spacetimes whose 
dynamics is reducible to scalar modes). In particular, since most formal and theoretical 
results obtained for LTB models can be readily applied to Szekeres models [62], the 
extension of our results to these models is currently under elaboration. In fact, the 
dynamical effect of q-scalars in LTB models may be applicable to even more general 
spacetimes. Specifically, we propose identifying fluid flow scalars that comply with 
FLRW dynamics (hopefully related to some weighed average), and then expressing 
local 1+3 scalars as fluctuations of the new scalars. The following step would be to 
explore if the 1+3 evolution equations given in terms of the new set of variables has 
the structure of evolution equations and constraints in the context of a formalism 
of exact perturbations on a FLRW background. We feel that this proposal is worth 
considering in future research. 



Appendix A. Analytic solutions. 

The metric function R in ([I]) can be obtained from the following parametric solutions 
of (2a) which define parabolic, hyperbolic and elliptic models: 



Parabolic: E = 0, R= (2Af) 1/3 rj 2 , t - t bb = - rf , (A.l) 

3 

M M 
Hyperbolic: E > 0, R = — (cosh r] - 1) , t-t bb = — ^ (sinh rj - ry) , (A.2) 

M M 
Elliptic: - 1< E < 0, R = (1 - cos r}) , t-t bb = (ry - sin 77) , (A.3) 

where t bb = t bb (r), which emerges as an "integration constant", is the "big bang time" 
because it marks the coordinate locus of the central expanding curvature singularity 



R(t,r) = for r 7^ constant. The classification of (A.l H A. 3) follows from the sign of 



E and g over ns the kinematic evolution of dust layers through the existence of a zero 



of R in (2a |: parabolic and hyperbolic layers perpetually expand and elliptic layers 
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expand/collapse. Since E = E(r), the models can exhibit, either the same kinematic 
pattern for all r, or "mixed" patterns that define models with regions of different 
kinematic behavior as E changes sign in specific ranges of r (see examples in 

The solutions ( A.l )-( A. 3 1 allow for a complete description of the dynamics of the 
models once we select the free functions M, E, t bb , though only two of these functions 
are strictly necessary because the metric ([I]) in invariant under re-scalings r = r(f) 
and thus any one of these free functions can be eliminated by a suitable choice of the 
radial coordinate. In particular, the expressions p6| |, Q, ([5| and §6§ for the basic 
covariant scalars p,W, fC, S, £ in |7| inv olve terms such as R! and R! , which can be 
computed by implicit derivation of (A.1)-(A.3). 



Appendix B. Functions of q scalars are q scalars. 



The relations ( 23 ) allow us to express local scalars A in terms of q-scalars A q through 



the differential inverse of the integral definition (13 1: 

(A)' g 



A = (A) q + 



V'/V q A « 



A' 



V'/V q 



(B.l) 



■ ql 'q 

where V' q /V q = 3R'/R and V q is defined in (94). As a consequence, any scalar 
expressible as a function of q-scalars is itself a q-scalar: if A q and B q are two q- 
scalars and U = U(A q , B q ), then 

3R' 
~R 



U' 



dA-y- A « ) + dB- {B - B « ] 



(B.2) 



x q u±j q 

so that we can identify Z q = U(A q ,B q ) as the q-scalar whose corresponding "local" 
scalar Z is given by 



(B.3) 



3R'/R 

where it is important to remark that (in general) we have Z =/= U(A,B). It is 
straightforward to show that Z and Z q above satisfy the integral definition ( 13 1, which 
basically states that a function of q-averages U((A) q , (B) q ) is itself the q-average (Z) q 



of a scalar Z given by (B.3 1, which in general does not coincide with U(A, B) nor its 
q-average (U(A,B)) q . 



Appendix C. Proof of (41) and (42) 



The proof of (41 1 for A = H is equivalent to the proof that the identity 



(Wji(r,ro))q[rb\ — holds for every domain D[rt,], where the scalar W% is given 
by: 

W n (r, r b ) = [D NL (H)] 2 - [D(H)] 2 - (H(r) - {H) q [r b ] f - (H(r) - H q (r) f . 



(C.l) 



-{'H)l[rb] holds, leading to 



We expand W n (r,r b ) and apply (21 1, so that (-ZH(r)(H) q [r b ] + {H)i[r b ]) q [r b ] 



(W n (r,r b )) q [r b ] = -(H) 2 q [r b ] + ^ {mn q - U\) V q df 



(C.2) 
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where V q is defined by (|94j), so that V' q = AttR 2 R' . Considering that (|3j) and Jl9 
imply 3% = 9 = V' q /V' q and 3V. q = 6 q = V q /V q . Inserting these identities into (C.2 
yields the desired result: 

1 r rb ( V 2 \ ' 

(W n (r,r b )) q [r b ] = - (H) 2 [r b } + / [^] df = -{U)\[r h ] + H 2 q (r b ) = 0, 



Vq(n) Jo 



(C.3) 



since {T-L) q [r b \ coincides with % q at the domain boundary r = r b . 



For A = p the proof of (41) is equivalent to proving W p (r, r b )) q [r b ] = 0, where 

\2 



W p (r, r b ) = [D NL (p)] 2 - [D(p)] 2 = (p(r) (p) q [r b ] ) 2 - ( p(r) p g (r) f . (C.4) 
Expanding W p (r,r b ), applying (21) and substituting p — M' /V' q and p q = M/V q from 



(2b) and (19), yields a relation similar to (C.2): 



1 f Tb 

(W p (r, r b )) q [r b ] = - (p) 2 q [r b \ + — — / (2pp q - p\) V q df 



(p) 2 q [n] 



M 2 

V 



dr 



-(p) 2 q [n} + p 2 q (n) = 0. 



(C.5) 



The proof of (42 ) for A, B = p, H follows along similar lines: it is equivalent to proving 
that: ( W p «(r, r b ) ) q [r b ] = holds, where: 

W pn (r,n) = D NL (H) D NL (p) - T>(H) B(H) 

= (P(r) - (P) q [n] )(H(r) - (U) q [r b ] ) - (p(r) - p q (r) ) (H(r) - H q (r) ) . (C.6) 
Substitution of the forms for "H, p, p q given before yields the desired result: 



(W P H(r,r b )) q [r b }=-(p) q [r b }(U) q [r ] 



= -(p) q [nW) g [r b ] + 



-— / {pH q + p q n - p q n q )V' q dr 

\ r b) JO 



V q (n) Jo 



MVg 

3V„ 



dr 



-(p) q [n}(H) q [r b }+ p q (r b )H q (r b ) = 0. 



(C.7) 



Appendix D. Proper curvature and kinematic tensors. 

We show in this appendix that all relevant (curvature and kinematic) proper tensors 
characteristic of LTB models are completely determined by 4 invariant scalars: the 
Ricci scalar, 1Z = 1Z C C , the only nonzero Newman-Penrose conformal invariant, \I/2, the 
Hubble expansion scalar W = 9/3 and the eigenvalue of the shear tensor, E. 
The Riemann and Weyl tensors are given by the following expressions: 



ab 



IK -cd — ^cd 



-U8^% + 28^11% 



C cd - d [c h d] 



2m r E „ u 



where square bracketed indices indicate anti-simetrization. The Ricci and electric 
Weyl tensors, TZ h d and E d follow from contraction of the tensors in (D.l ): 

K b d = K b cd = \ll (h b d + u b u d ) , 



E\ = u a u c Cf c 



(D.2) 
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where e b d = 8 b d — 3n b n d is the common eigcnframc for E h d and the shear tensor a d 



defined in ^ and (|6|. By combining (D.ll and (D.2| we get 



c u d] 

Cf d = * 2 (*j-ej + 2u [c e^) = * 2 (/£ 3u [c u^) e% 



°[c°d] 



(D.3) 
(D.4) 



Hence, all curvature proper tensors are determined by two invariant scalars: 1Z 
and \?2) which relate to the density and its local fluctuation by 1Z — 8irp and 
*2 = -£ = (4tt/3)(p- Pq ) = (4tt/3)D(p). 

The kinematic proper tensors are the shear and "expansion" tensors: o~ a b, Hat- 
Both are expressible in terms of the invariant scalars H = 8/3 = V Q it a /3 and 



o-abe 



ah . 



0~ab 



Hab — 



(D.5) 



where the explicit forms of ri and S are given in ^ and ^. Notice that £ = -T>(U) 
is a local fluctuation of ri. 
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